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Abstract 

In this paper we study a renormalization-group map: the block averaging 
transformation applied to Gibbs measures relative to a class of finite range lattice 
gases, when suitable strong mixing conditions are satisfied. Using block decimation 
procedure, cluster expansion (like in [HK]) and detailed comparison between statis- 
tical ensembles, we are able to prove Gibbsianess and convergence to a trivial (i.e. 
Gaussian and product) fixed point. Our results apply to 2D standard Ising model 
at any temperature above the critical one and arbitrary magnetic field. 
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1. Introduction 

This paper concerns the rigorous analysis of some real-space renormalization group trans- 
formations (RGT) (see, for instance, ref. [NL] for a general introduction to this subject). In the 
recent years many works have been devoted to the question of well-definedness of RGT. We refer 
to ref. [EFS] for a clear and complete discussion of the problematic as well as for an exhaustive 
description of the general setup of renormalization maps from the point of view of rigorous sta- 
tistical mechanics. Already in the seventies (see [GaK], [CG], [GP], [I]) the question was raised 
of whether or not some typical RGT give rise to a well defined renormalized interaction. In other 
words, calling fi^ the renormalized measure arising from the application of a RGT on scale i to 
the Gibbs measure y, we pose the question of Gibbsianess of fi^ , namely we ask ourselves whether 
is a Gibbs measure corresponding to a finite-norm translationally invariant potential so that 
the "renormalized Hamiltonian" is well defined. 

More explicitly: let us assume that our RGT can be expressed as 

^V) = £^v>)/^) (i.i) 



where T^(o~',o~) is a normalized non-negative kernel. The system described in terms of the a 
variables by the original measure fi is called object system. The cr"s are the renormalized variables 
of the image system described by the renormalized measure fj,^ . We can think of the transformation 
Ti as directly acting at infinite volume or we can consider a finite volume version and subsequently 
try to perform the thermodynamic limit (see [EFS] ). 

The above mentioned pathological behavior (non-Gibbsianess of //^) can be a consequence 
of the violation of a necessary condition for Gibbsianess called quasi-locality (see [Ko] , [EFS] ) . The 
latter is a continuity property of the finite volume conditional probabilities of fi^ which, roughly 
speaking, says that they are almost independent of very far away conditioning spins. In many 
interesting examples (see [El], [E2], [EFS], [EFK]) violation of quasi-locality and consequently non- 
Gibbsianess of the renormalized measure fj,^ is a direct consequence of the appearance of a first 
order phase transition for the original (object) system conditioned to some particular configuration 
of the image system. More precisely, given a configuration a', let us consider the probability 
measure on the original spin variables given by 



It defines the constrained model corresponding to a' (which here plays the role of an external 
parameter). For some particular a' it may happen that the corresponding measure /i CT '(<r) exhibits 
long range order inducing violation of quasi-locality and then non-Gibbsianess for the image system. 
See [EFS] and also [GP], [I] where this mechanism was first pointed out. 

This pathological behavior is often induced by configurations a' highly non-typical with 
respect to the measure yS e \ This suggests the introduction of a weaker notion of Gibbsianess 
requiring well-definedness of renormalized interactions not for all renormalized configurations a' 
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but, rather, for ^ -almost all a 1 (see [D2], [BKL], [LM]). However also this point of view poses 
various other problems (see [ES] and references therein). 

It is also natural to ask ourselves about robustness of this pathology. Sometimes it can be 
shown that non-Gibbsianess is an artifact due to a wrong choice of the scale £ of the transformation 
in terms of the thermodynamic parameters of the object system. For instance in [M04] has been 
considered the case of the measure fj,^ = T^p lia,hi where T^p is the so-called decimation transfor- 
mation on scale £ (see [EFS]) and \ip,h is the Gibbs measure for the standard 2D Ising model, h and 
(3 being, respectively, the external field and the inverse temperature. In [EFS] the authors show 
non-Gibbsianess for some choices of h,(3,£. On the contrary, in [M04] it is shown that, for the 
same values of h, [3 for which, for suitable £, in [EFS] the authors got non-Gibbsianess, by changing 
£ into a sufficiently large £' , depending on (3,h, one gets again Gibbsianess. We also mention it 
is possible to show that, by iterating the transformation, one has convergence to a (trivial) fixed 
point, see [M04] and also [Ka] for the high temperature case. The above behavior is related to the 
fact that, given suitable values of the parameters (3, h (close to the coexistence line h = 0,(3 > (3 C ), 
on a suitable scale £ some constrained models can undergo a phase transition (somehow related 
to the phase transition of the object system); whereas, given the same h,j3, for sufficiently large 
scale £ any constrained model is in the weak coupling region. Another notion of robustness of the 
pathology is related to the application of decimation transformations, see [LV], [M05]. 

Let us stress that the fact that the object system is very well behaved in the sense that it is 
in the unique phase region (in the strongest possible sense) does not preclude the possibility that 
some constrained model undergoes a dangerous phase transition inducing the pathology. 

On the positive side, since the pioneering paper [CG], there are many indications that if the 
constrained models are in the weak coupling regime, then Gibbsianess of the renormalized measure 
follows. Recently Haller and Kennedy gave very interesting new rigorous results in this direction. 
They proved, under very general hypotheses, that if all constrained models satisfy a uniform (in the 
constraint) version of the Dobrushin-Shlosman complete analyticity condition (see [DS2], [DS3]) 
then the renormalized measure is Gibbsian with a finite norm potential which can be computed via 
a convergent cluster expansion. 

Another interesting question, which, in a sense, is the main object of the present paper, is 
the convergence of the iterates of RGT or, in other words, the behavior of the transformation 
for large I. This problem has not been, up to now, studied very much from a point of view of 
rigorous statistical mechanics. Here we present results referring to non-critical systems and so we 
have convergence to a trivial fixed point, i.e. Gaussian and product (which correspond to infinite 
temperature). Indeed most of the recent papers concerning rigorous results on RGT refer to the 
non-critical region with some exceptions, see [BMO], [CiO], [HK], where the authors consider 2D 
critical Ising system but only for one step of RGT. 

Let us now introduce, for the standard 2D Ising model, the Block-Averaging Transformation 
(BAT). It is convenient to use the lattice gas variables. For a standard Ising system enclosed in a 
finite volume AcZ 2 the configuration space is therefore {0, 1} A ; given a configuration 77 G {0, 1} A 
and a site x G A, r) x G {0, 1} represents the occupation number at x. For free or periodic boundary 
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conditions the energy associated to a configuration n £ {0, 1} A is: 

E A ( V ):=-(3 Wy-X^V* (1-2) 

(x,y)CA x€A 

where (x, y) is a pair of nearest neighbor sites, (3 is the inverse temperature and A is (3 times the 
chemical potential so that the Boltzmann factor is ~ exp{— £^a (??)}• Given (3 let A* = X*((3) be 
the value of A corresponding to the value zero of the magnetic field h appearing in the expression 
of E in terms of spin variables a x = 2n x — 1. For (3, A in the uniqueness region: ((3, A) € {/? < 
/3 C } U {/3 > /3 C ,A / A*} (/3 C is the inverse critical temperature), let p,p t \ be the unique infinite 
volume Gibbs measure. We partition 1? into square blocks Qi(i) of side i and centers at the 
points i belonging to the rescaled lattice (^Z) 2 . Let iV, = Ni(w) := YlxeQeii)^ ^ e ^ ne num b er 
of particles in the block Qg(i) in the r\ configuration, p = p((3,X) = p/3 t \(r] ) be the equilibrium 
density, x = x(A A) ■= J2 x€ z 2 [PpAvoVx) ~ Vf3,x(rio)Vf3,\{Vx)] be the susceptibility; we then set: 

Mi := N '-pW<\ (1.3) 
the random variables Mj are centered and normalized; they take values in 

n(0 ._ ( -p\Qi\ i-p\Qi\ \Q t \(i-p) \ (14) 
1 ' IvIWvlto''"' VWx J 

We expect Mi, i e (£Z) 2 to have a product (Gaussian) limiting distribution as i — > oo. 

The renormalized measure = /U^ A (arising from the application of the BAT transforma- 
tion on scale I to /U^a) is the joint distribution of the random variables Mj's under yU/^Aj i- e - it 
is obtained by assigning to each block Qg(i) a value rrii £ and by computing the probability, 
w.r.t. the original Gibbs measure pp^\ of the event: Mj(r/) = rrii. In other words, in the notation 
of (1.1) in the case of BAT we have: a = {n x },a' = {rrii} and 



T ( B e) AT (m, V ) 



1 if Mi(r]) = m, Vi 
otherwise 



In this case a constrained model is a multi-canonical Ising model; namely an Ising model subject 
to the constraint of having a fixed number of particles in each block Qg(i). 

Theorem 1.1. Consider a 2D Ising system with [3 < (3 C and A € R given. Then there exists 
£o 6 N such that Ml > lo is Gibbsian with a finite norm translationally invariant potential 

$W = {$^ } , X C («) d }. 

Furthermore it is possible to decompose the potential into a short and a long range part, 
$M = $(/),«• + $(/),ir j w]iere 3k g n . = o if diam(A) > « and we have the following: 

fi) there is a > such that 



lim V e Q l x l sup $; ),lr (m„ia) 



X90 m,6Si; 
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(ii) there exist a > such that 



lim sup ^^' sr (rrii,i & X) 

l— >00 - (I) 



lim sup 



{/} ( m i) + 



for \X\ > 2 



for i £ 



We want to stress that the results hold for £ sufficiently large. Certainly, in particular, we 
cannot exclude that, very near to T c , for some, not sufficiently large £, the renormalized measure 
is not Gibbsian. On the other side, it is easily seen that taking the limit £ — > oo is equivalent to 
iterate BAT transformation on a given scale £o; to show this it is sufficient to take £ = £q with 
n G N. Theorem 1.1 above says that not only the renormalized measure /A \, for any sufficiently 
large £ is Gibbsian but the corresponding renormalized potential actually converges, as £ — > oo, 
to the one of a system of non-interacting harmonic oscillators. 

We notice that the limiting image system as £ — > oo becomes an unbounded spin system and 
the usual setup of Gibbsianess does not apply to it (see [EFS]). It is therefore clear that we have to 
introduce a large field cutoff. Indeed our result is almost optimal as we introduce this cutoff only 
for the short range part of the interaction and, moreover, the cutoff diverges as a power law in £. 
On the other side it is not difficult to convince ourselves that the convergence result, at least in the 
form given above, cannot hold without any restriction on the large fields. 

This paper contains also other, much weaker, results that apply to Ising model below T c at 
A A*, see Theorem 2.3 below. In that case we are forced to restrict the possible values of m; also 
in the computation of long-range part of the renormalized potential; indeed we have of course to 
forbid that rrii lies in the phase coexistence interval. 

Results in the same direction as Theorem 1.1 were obtained by Cammarota [C]; the main 
differences w.r.t. the present paper are that Cammarota considers a high temperature (much higher 
than T c = f3~ v ) situation and that he introduces a finite (not growing to infinity as £ — > oo) field 
cutoff. The approach of [C] is substantially different w.r.t. ours; [C] uses a high temperature 
expansion: the small parameter is (5 and the system is supposed to be weakly coupled on scale 
one; whereas since we want to treat a system with T = j3~ l higher but arbitrarily close to T c , 
we have to use an approach supposing weak coupling only on a sufficiently large scale depending 
on the temperature T > T c that we have chosen; indeed we are forced to act in the scenario of 
the so called restricted complete analyticity. Let us try to clarify this point. Exactly in the spirit 
of renormalization group theory we can say that a system above its critical point is very weakly 
coupled on a scale large compared to the correlation length; as we want to consider any T > T c 
we have to take into account the divergence of the correlation length when approaching T c (from 
above). The above statement: "the system is weakly coupled on a scale larger than the correlation 
length" seems a tautology; in fact it is not since we need a suitable mathematical setup in order 
to be able to implement the above simple observation. The basic idea is to obtain a perturbative 
expansion on the basis of very strong mixing conditions satisfied by the Gibbs measures; the small 
parameter ceases to be the inverse temperature but it will, rather, be related to the ratio between 
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the correlation length and the scale on which we are analyzing our system. The geometrical objects 
(polymers) in terms of which we perform our perturbative expansions will not live any more on 
scale one (like in [C]) but on a scale sufficiently larger than the correlation length. 

A possible notion of strong mixing is the exponential decay of truncated correlations for any 
finite volume Gibbs measure with decay constants uniform in the volume and in the boundary 
conditions. This is a stronger notion w.r.t uniqueness of the Gibbs state or exponential clustering 
of infinite volume truncated correlations. In [DS2], [DS3] Dobrushin and Shlosman introduced and 
studied the so called completely analytical interactions showing, in those cases, the above strong 
mixing behavior for any finite or infinite domain of arbitrary shape. This complete analyticity turns 
out to be a too strong notion in the context of renormalization group theory. Indeed Dobrushin- 
Shlosman's complete analyticity implies that exponential clustering takes place even inside volumes 
with very anomalous shapes (for instance with anomalous ratio between boundary and bulk) so that 
one is forced to take into account the influence of boundary conditions up to a scale of order one. 
There are cases of systems perfectly well behaved on regular domains, say cubes, which, however, 
do not satisfy D-S complete analyticity because of their behavior for anomalous shape domains (see 
example in [M02]). Another point of view, introduced in [O], [OP], [MOl], [M02], [M03] leading 
to what can be called restricted complete analyticity, takes into account only regular domains. In 
this approach there is a minimal basic length L and one never goes below L in the sense that one 
only considers domains obtained as disjoint unions of cubes of side L (for instance cubes of side 
nL). 

The algorithm used in the present paper to compute the renormalized potential is the fol- 
lowing. We start, as basic hypothesis, from restricted complete analyticity for the constrained, 
multi-canonical systems, with a minimal length L proportional to the scale i of our BAT transfor- 
mation and with decay constants uniform in the constraint. We then construct a convergent cluster 
expansion which allows us to compute the renormalized potential. Since studying directly the mix- 
ing properties of a canonical or multi-canonical measure is a very difficult task we instead deduce 
it by using a sharp form of equivalence, or better comparison, between canonical and grancanonical 
ensembles. Indeed, the main key novel technical point of this paper is to get a very precise notion 
of equivalence of ensembles, implying the validity of a finite size condition, which, in turn, will 
imply a strong mixing condition for the constrained multi-canonical systems. See also [DT], [CM], 
[Y] for a further discussion on the equivalence of ensembles. 

Certainly assuming strong mixing for the object system with a given value A of the chemical 
potential is not sufficient to imply the strong mixing property of the constrained models even 
at the level of regular domains. It is, rather, necessary to assume for the object grancanonical 
system a strong mixing condition uniform in A. Quite surprisingly, this condition is not sufficient 
in general. Indeed it turns out that what we really need is a strong mixing condition for a multi- 
grancanonical object system; by multi-grancanonical we mean a grancanonical measure which is not 
translationally invariant because in each cube Qe(i) we put a different chemical potential A^ whereas 
we leave the original, translationally invariant, mutual interaction. It happens, as it is shown by 
an example in Appendix A. 2 that uniform (in A) strong mixing for a grancanonical measure does 
not imply uniform in A = {A^} strong mixing for the multi-grancanonical measure; this pathology 
is due to the possibility of a sort of layering phase transition, with long range order, taking place 
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along the interface between two contiguous large cubes with different chemical potentials Ai,A 2 
even though, introducing the same chemical potential Ai = A2 = A in both cubes the resulting 
system, is, V A, very well behaved. On the other side we show that since the interface between two 
regular two-dimensional domains is one-dimensional, this layering phase transition cannot occur 
when the object system lives in two-dimensions. Then the result of Theorem 1.1 ultimately follows 
from strong mixing, uniform in A, exploiting the two-dimensionality of the Ising system. The latter 
follows from the general result of [MOS] saying that in two dimensions the so-called weak mixing 
implies strong mixing, provided one is able to prove weak mixing for our particular model. This 
is a weaker notion of mixing of a finite volume Gibbs measure saying, roughly speaking, that the 
influence of a change in a conditioning spin on a site x outside a domain A decays, inside A, with 
the distance from the boundary dA and not, like it would be the case assuming strong mixing, 
with the distance from x. Weak mixing, uniform in the chemical potential A, for Ising model above 
T c has been proved by Higuchi in [H] exploiting general results by Aizenman et al. [ABF] about 
boundedness of susceptibility above T c . We thus use the two-dimensionality in two crucial points: 
i) in deducing uniform strong mixing for multi-grancanonical measure from the same property for 
simple grancanonical measure and ii) in deducing strong mixing from weak mixing. On the other 
side, given the strong mixing condition for the multi-grancanonical measure, the results on the 
RGT of this paper apply to any dimension. 

The general results about Gibbsianess of renormalized measures that have been obtained by 
Haller and Kennedy in [HK], use a strategy very similar to ours. Indeed their computations, also 
based on the methods developed in [O] , [OP] are much simpler and more transparent then ours but 
apply only to the case when the image system is Isingdike; namely the a' variables are dichotomic 
as in Majority rule or Kadanoff (see [EFS]) transformations. Haller and Kennedy for a given £ use 
the hypothesis of D-S complete analyticity of the constrained models to deduce Gibbsianess of the 
measure resulting from the application of one transformation. 

We conclude by a brief outline of the various steps needed in the proof of Theorem 1.1. Higuchi 
[H] proves weak mixing, uniform in A, for Ising model above T c . [MOS] proves, in general, that in 
two dimensions, for regular domains, weak mixing implies strong mixing. In Appendix A.l we prove 
that in two dimension strong mixing uniform in A for grancanonical measure implies strong mixing 
uniform in X = {Ai} for the corresponding multi-grancanonical measure. In Section 4 we prove 
results about comparison, in a finite volume A, between multi-grancanonical and multi-canonical 
measures with precise estimates of the behavior in A. From this and previous points we deduce 
that, on a sufficiently large scale, an effective (propagating to arbitrarily large, regular domains) 
finite size condition is satisfied for multi-canonical constrained systems. Then, from this finite 
size conditions, using the theory developed in [O], [OP] we are able to perturbatively compute the 
renormalized Hamiltonian; and to extract the potentials. The long range terms of the interaction 
potential are computed starting from a cluster expansion whose convergence is directly related to 
the validity of the above finite size condition. Finally the short range terms are handled via a local 
central limit theorem for the multi-grancanonical measure. 
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2. Notation and results 

We introduce the general setup: the one of the finite state space, lattice spin systems. Contrary 
to the usual treatments we drop the hypothesis of translation invariance; indeed it will be replaced 
by spatial uniformity of some basic estimates. We start by giving a list of basic definitions. 

2.1. The lattice. 

For x = (x\, • • • ,Xd) G K d we let |x| := sup fe=1 ... jd \xk\- The spatial structure is modeled by the 
d-dimensional lattice C := Z d in which we let ei,i = 1, . . . ,d be the coordinate unit vectors. We 
shall denote by x, y, ■ ■ ■ the sites in C and by A, A, • • • subsets of C. We consider C endowed with 
the distance d(x, y) = \x — y\. We use A c := C \ A to denote the complement of A. For A a finite 
subset of C (we use A CC C to indicate that A is finite), |A| denotes the cardinality of A. For 
x G C and I an odd integer we let Qe(x) := {y G C : d(y,x) < (i — l)/2} be the cube of side £ 
centered at x; for £ an even integer we let instead Qe(x) := {y G C : \y — (x + e)| < 1/2}, where 
e := (1/2, • • • , 1/2), be the cube of side £ centered inx + e (which belongs to the dual lattice). We 
shall denote Qe(0) simply by Qi. Given r > and A C C we introduce the outer boundary of A by 
d r A := {x G" A : d(x, A) < r}. We let also A" := AU d r A. 

Given an integer £, we also introduce the rescaled lattice Cg := (£Z) d which is naturally 
embedded in C; we shall therefore regard points in also as points in C without further mention, 
more precisely we will make the following identification: Cg 3 . . . , id) = ■ ■ ■ , £id) G £>■ We 
use i,j, ■ ■ ■ to denote points in C( and /, • • • to denote subsets of Cg. Analogously the distance in 
£i is denoted by dg(i,j), therefore for i,j G Cg we have d(i,j) = £dg(i,j). 

2.2. The configuration space. 

We suppose given a positive integer N G N + and, for every x G £, a positive integer N x < N. We 
then introduce the following: 

- Configuration space of a single spin. For any x G C we have a finite set £l x , \Q X \ = J\f x + 1. 
We identify Vt x with {0, 1, . . . ,J\T X } which we consider endowed with the discrete topology; 

- Configuration space in a subset A C C. We set '■= CSxeA^x! 

- Configuration space in the whole C. We set := ® x ^c^x and equip it with the product 
topology. 

We can therefore look at a configuration a G as a function a : C i— > {0, 1, . . . ,J\f}. The integer 
a x = a(x) is called value of the spin at the site x G A in the configuration a. For A C C, we use 
oa := {a x G fl x , x G A} to denote the collection of spins in A. For x G C we define the shift -Q x 
(acting on by {$ x a)y := cj^-^. 

We also introduce C (fi^) the space of continuous functions on f2 which becomes a Banach 
space under the norm ||/|| := sup CT |/(c)| and note that the local functions (i.e. the functions 
depending only on a finite number of spins) are dense in C (fi^)) . For / a local function depending 
on the spins in A CC £, i.e. f(a) = /(<ta), we let S(f) = supp(/) := A be the support of /. 

In the case Af = 1 the spin a x takes values in {0, 1}, i.e. we have a lattice gas. In such a case 
we use the notation O := {0, 1} £ and denote by 77, (, ■ ■ ■ typical elements of fi; the value r\ x G {0, 1} 
is interpreted as the occupation number in x. Given r\ G O we define a new configuration n x which 
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is obtained from rj by flipping the occupation number in x, i.e. 



(Vly ■■- 



Vy \iy + x 
l-r] x if y = x 



2.3. The potential. 

A potential <I> = {<£a,A CC £} is a family of functions labeled by finite subsets of C and 
<I>A : ^a^ ^ ^- We introduce the following possible conditions on <E>: 

- Finite range. There exists r > such that <J>a = if diam(A) > r; 

- Translation invariance. For each x 6 C, $a(<7) = &A+x($xO~)- 

We note that the potentials (which do not need to satisfy the conditions above) form a linear 
space. Given a > 0, we introduce in it the norm || • || Q defined by 

||$||c := sup Ve a|A| sup |$a(ct a )| 

We also note that in the translation invariant case we can omit the first supremum above. 

Given A CC £ and a potential <1? with bounded || • 1 1 norm, the energy associated to a 
configuration a when the boundary condition outside A is (the restriction to A c of) r G SI, is given 
by: 

E a {*\t) := $r(<T° A r) (2- 1 ) 

rnA/0 

where 

a x if x £ A 



(a o A r) 



t x if x A 



Note that the sum on the r.h.s. of (2.1) is absolutely convergent (uniformly in a and r) by the 
boundedness of ||3>||o- We also remark that for a finite range potential the map r i— > E\(o\t) 
depends only on TQ r A- 

2.4- The Gibbs measures. 

Given a potential $ of bounded || • ||o norm, for each A CC £ we define the (finite volume) Gibbs 
measure in A with boundary condition r as 

Ma(^) : = ^t^v{-Ek{(t\t)} 

Z A 

where Z^, called partition function, is the normalization constant, i.e. 



Zl = Zl{^):= Y, e 



Ea(ct\t) 

(AT) 



Note that we have included the inverse temperature in the definition of energy; in fact it will be 
kept fixed in our analysis. 
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We regard \x T K also as a measure on the whole by giving zero mass to the configurations 
a which do not agree with r on A c . The (infinite volume) Gibbs states associated to the potential 
$ are then the measures [i on which satisfy the DLR equations 

J /i(dr) M X(/) = /*(/), for any A CC C, f € C (ft (A °) 

For a translationally invariant lattice gas we observe that we have <&{ x }(rj) = —Xr] x + a for 
some constants A, a € R. We neglect the constant a (which do not affect the definition of the Gibbs 
measure) and note that A is interpreted as the chemical potential. We also introduce the activity 
z e R + by z := e A which we use to parametrize lattice gases with different chemical potentials. In 
such a case we write $ = (z, U) where U = {^a, A CC C, |A| > 1} and call U the interaction. We 
shall also write (sometimes) fx\ z (resp. Z\(z)) in order to indicate explicitly the dependence on 
the activity z. 

2.5. Strong Mixing Conditions. 

In what follows we recall notions concerning some mixing properties of Gibbs measures. Most of the 
theory has been, up to now, developed in the finite range, translationally invariant case. Extension 
to not translationally invariant cases, when suitable uniform conditions hold, is, in most of the cases, 
straightforward. In particular we will be concerned with the so-called strong mixing condition which 
can be formulated in terms of exponential clustering of truncated expectation with respect to the 
Gibbs measures in certain domains A with r boundary conditions when this exponential clustering 
takes place uniformly in A and r. This strong mixing condition implies uniqueness of infinite volume 
Gibbs measure and its exponential clustering. It can be shown that finite volume strong mixing 
condition, with constants uniform in the volume and in the boundary conditions, is strictly stronger 
than the equivalent infinite volume notion (see [Sh], [Ba], [DM], [CM]). As it has been shown by 
Dobrushin and Shlosman (see [DS2], [DS3]) this strong mixing condition, supposed to hold for any 
(finite or infinite) volume A, is equivalent to many other conditions like analyticity properties of 
thermodynamic and correlation functions or tree-decay of semi-invariants. Interactions giving rise 
to this kind of nice behavior have been called by Dobrushin and Shlosman completely analytical. 
Among their equivalent complete analyticity conditions, Dobrushin and Shlosman have introduced 
suitable finite size conditions that they call "constructive conditions". They show that, supposing 
that there exists a finite domain A such that strong mixing condition is satisfied with suitable 
(depending on A) decay constants for all subsets of A, then a strong mixing condition holds for all 
(finite or infinite) volumes. 

We refer to [M02] for a discussion on the applicability of this point of view. Indeed often 
the request of exponential clustering for arbitrary shape does not fit with many reasonable appli- 
cations. There are examples (see [MOl]) where nice exponential mixing properties hold for regular 
domains (like, for instance, cubes) and in infinite volume, whereas they fail to hold for domains with 
anomalous ratio between boundary and bulk, implying violation of Dobrushin-Shlosman complete 
analyticity. In [O], [OP], [MOl], [M02], [M03] another scenario has been introduced, more suited 
to the renormalization group problematic. It can be called "restricted complete analyticity" or 
"complete analyticity for regular domains" This point of view refers to exponential mixing in finite 
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volumes multiples of a given cube of size Iq. In the framework of this theory one can develop a 
constructive condition of the following kind: if 3£q such that a suitable (depending on Iq) mixing 
condition holds in the cube Qi , then the same condition (possibly with worse constants) holds 
for any multiple of Qe . This possibility of propagation from finite to arbitrarily large (and even 
infinite) volumes is called "effectiveness" in [M02]. Subsequently many results have been obtained 
in the framework of restricted complete analyticity that could have been problematic and even 
false in the context of Dobrushin-Shlosman complete analyticity (see, for instance [M02], [M03], 
[MOS], [SS]). 

Given a measure [i and two square integrable random variables f,g we denote by fi(f;g) ■= 
~ M/VG?) the- covariance between / and g. For A CC £ we introduce ^X ; A as the rela- 

tivization (projection) of the Gibbs measure /x A to i.e. 



^a ; a(o"a) := J lAidC,) 



L Ca=o"a 



We finally recall that the total variation distance between two measures /x, v on a finite set S is 
given by 

Var( M) v) := \ £ |mM ~ = sup | M (X) - v(X)\ 

If a, b € M, we let a A 6 := min{a, b}. For a finite range potential we introduce the following strong 
mixing condition. 

Condition SM(£q) (Strong Mixing). 

Given an integer £ we say that the potential $ satisfies SM(£ ) if there exist two constants A, 7 > 
such that for any volume 

A=(jQe (i), ICC£ eo (2.2) 
iei 

the following bound holds. For any x € d r A and any A C A we have 

sup sup Var f/x A ° A a , ^. A ) < Ae~^ x '^ (2.3) 



We next discuss finite size conditions which imply SM(^o)- Let m be an integer, m > r, and 
consider the cube Q3 m (j), j £ C m . Given a particular lattice direction a we can partition Qz m (j) 
into three parallelepipeds having d — 1 sides equal to 3m and the last one equal to m (slices) with 
the minimal side parallel to the direction (slice orthogonal to e^). We write 

Qsmti) = Qt~\j) u QL 0) (J) u (2-4) 

here Q^\j) denotes the central slice. 

Let Pm(j) be the set of all subsets of Q^\j) which are unions of cubes Q m (j)- For V € 
Pm(j) let d^'^V, d^'~^V denote the part of d r V in the direction of e,, — e, respectively (opposite 
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r-faces of V). Given a, (,t G ft, y G d^'^V, y' G d^' ^V, we denote by a^' + \C > ^' \t the 
configuration obtained from r by substituting in d^'^V fl^'^V the restrictions of a, C, respectively: 



' a x iixed^V 
t t otherwise 



analogously we denote by a y , C* 1 ' \t the configuration obtained from r by substituting to r in y, 
d^'~^V the restrictions of a,(, respectively: 

( &x i£ x = y 

I t x otherwise 

finally we denote by a y , ( y >,r the configuration obtained by substituting to r in y, y' the restrictions 
of a, (: 

( a x \ix = y 

K>Cj/,7-) x == (x if x = y' 



of course 



^ otherwise 



T (i, + ) )T (t, ) >T = TwjT (*. \ T = Ty,Ty>,T = T. 



We introduce the notation: Zy(r) := Z v . 
Condition Cl(m,£i) (see /OP/, Pg fP<^- 



sup sup sup sup 
j€C m i€{i, ■■■,<!} veP^\j) °-,TenW) 

Condition C2(m, e 2 ) 



Z y (^+),q^-),r) Z y (r^+),r^-),r) 
Z v (o-( i '+),r( i -),r) Z y (t( < >+),(t( < -),r) 



< El 



(2.5) 



sup sup sup sup sup 
je£ m ie{i, ■■■,<!} veP^U) y^^+W <7,TefiM 



(o-j,, r(* •"), t)Z v (r y , cr(» >"), r) 



< 



£2 



I (2-6) 



Condition C3(m, e 3 ) (see [OJ). 



sup sup sup sup sup 

jeC m ie{l,—,d] VGP^d) yedi^+lv ff.ref!^) 



Z v (a y ,a y >,T)Z v (T y ,T y/ ,T) 



Zv(<r y ,T y >,T)Z v (T y ,a y ,,T) 



< 



£3 



m 



2(d-l) 



(2.7) 



It is easy to show, using a telescopic argument, that there exists a constants k such that 
C2(m,e2) implies Cl(m, k£ 2 ) and C3(m, £3) implies C2(m, k£ 3 ) (see [O]). It is also immediate to 
see that the results proven for the translationally invariant case in [O], [OP], [M02] extend to the 
general case when the space uniform condition holds. We have indeed the following result. Let 



e(d) := [3(2 d+1 + l)]" d 2~' Zd e" 4 , 



d r,-2d -4 



(2.8) 
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then condition Cl(m, e(d)) implies the existence of a convergent cluster expansion which, in turn, 
implies SM(m). 

We remark that once we have proven the crucial point which is the effectiveness, namely that 
Cl(m, e(d)) implies SM(m), then, considering the rescaled system whose new single spin variables, 
labeled by j £ C m , are the old spin configurations in the blocks Q m (j)i we can apply Dobrushin- 
Shlosman's results [DS1], [DS2], [DS3] to get all their equivalent mixing and analyticity properties 
of the Gibbs state for every "multiple" of the Q m 's namely for all volumes A of the form (2.2). 
This is the restricted complete analyticity namely the validity of the D-S equivalent properties (see 
[DS2], [DS3]) for every volume of the form (2.2). In particular SM(£ ) is equivalent to: 

Condition SM2(^ ) 

Given an integer £q we say that the potential <1? satisfies SM2(£q) if there exist two constants 
A, 7 > such that for every pair of local functions f, g and every volume of the form (2.2) 

sup \^ A (f;g)\<A(\S f \A\S g \)\\f\\ ll^ll e-^ s " s ») (2.9) 

where we recall Sf,S g are the supports of f,g. 

Indeed the implication Cl(lo,e(d)) SM2(^o) can be obtained directly via cluster expansion 
by using the methods of references [O], [OP]. We do not reproduce here the results of [O], [OP] 
but, looking at the application to the renormalization group problem that will be developed in next 
section, the reader could easily understand these results. 

2.6. Lattice gases: Uniform Strong Mixing Conditions. 

We here consider just a finite range lattice gas (i.e. £1 = {0, 1} C ) and introduce some uniform strong 
mixing conditions which are needed to study the RG map. These conditions say - rougly speaking - 
that SM(4>) holds uniformly in the external field (one body interaction). Unfortunately, as discussed 
in the Introduction, we need such a condition also for some non homogeneous external field. Such a 
condition plays also a crucial role in the ergodic properties of the Kawasaki (conservative) dynamics, 
[Y]. 

Given a finite range lattice gas with translationally invariant interaction we introduce the 
following Condition. We recall z = exp{A} is the activity. 

Condition USM(„4) (Uniform Strong Mixing). 

Given an open set A C [0, oo), we say that the interaction U satishes USM(A) if for each z G A 
there exists Iq = £q(z) such that condition SM(£ ) holds for (z, U). Furthermore the following is to 
be satished: 

(i) for any closed set C C A we can take the constants £q,A,^ uniform for z £ C; 

(ii) we can take A = Aqz A z _1 for some other constant Aq independent of z. 

Remark. We note that for A = [0, e] U [e _1 , oo) with e small enough (depending on d, r and ||3>||o) 
the above conditions holds. Indeed for z A z~ x small, SM(1) follows from standard perturbative 
theory (for instance by using Dobrushin single site criterion [Dl]). We can therefore safely replace 
the set [0, oo) in Condition USM(^l) by the compact set [e, e -1 ]. To avoid delicate continuity 
questions we introduced (i) above as an independent hypothesis. The same argument shows (ii) is 
automatically satisfied; we have included it only for convenience. 
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Condition GXJSNL(Global Uniform Strong Mixing) 

We say that Condition GUSM is satisfied if Condition USM(A) holds for A = [0, oo). 



2.7. The multi-grancanonical state. 

Let I be a positive integer and A C C a disjoint union of cubes of side £, i.e. 



A=\jQi(i), ICC £ (2.10) 

iei 

Given a lattice gas with a finite range translationally invariant interaction U, we next define a 
Gibbs measure in $7 a which has a fixed chemical potential in each cube Qe. We call such a measure 
a multi-grancanonical state. Let z := {zi £ [0, oo), i £ 1} the measure /x A z is then defined as a 
Gibbs measure in £l\ whose potential <£- = (z, U) is given by 

z_, . _ J -77s log ^ if T = {x} and x G 

rW) "~K(i) if|r|>i 

If I CC the finite volume multi-grancanonical measure is thus defined by 



iei 




^ {V) = Zl(zj • ex P 1 E $ r(r? o A r) } (2.11) 



|T|>1 



where Z T K (z) is the normalization constant and 

xeQe(i) 

is the total number of particles in Qe(i). We stress that the multi-grancanonical state ijl t k z does 
depend on I. 

We shall need a stronger version of Condition USM which is formulated in terms of the 
multi-grancanonical state. 

Condition MUSM(yl) ( Uniform Strong Mixing for Multigrancanonical States). 
Given an open set A C [0, oo), we say that the interaction U satisfies MUSM(A) if the following 
condition holds. For each closed set C ^ A there are constants £q 6N,i,7>0 such that for any 
t integer multiple of any / CC and any z G C 1 we have that for any A of the form given in 
(2.10) the multi-grancanonical measure \i\ z satisfies the bound (2.3). 

Condition GMUSM (Global Uniform Strong Mixing for Multigrancanonical States). 

If Condition MUSM(A) holds for A = [0, oo) we finally say that Condition GMUSM is satisfied. 

We also give an effective finite size condition of type CI which implies MUSM(^l). We note 
that if V € P$(j) we have V = U k( _yQ m (k) for some V C C m uniquely determined by V. We 
denote below by Zy ;£ (r) the multi-grancanonical partition function as defined in (2.11) with I = m. 
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Condition MUCl(i) 

Given an open set A C [0, oo] we say that MUC1(A) holds for the interaction U if for each closed 
set C C A there exists an integer m such that 




- 1 < e(d) (2.13) 



Indeed, by exploiting the translationally invariance of U and following the same argument as 
the one used in [O], [OP] it is easy to verify that if MUC1 holds we have that also MUSM(.A) holds. 

Remark 1. In the high temperature regime, ||£7||o < £ with e small enough, it is not difficult to 
show (by using, for instance Dobrushin's single site condition [Dl]) that Condition MUC1([0, oo)) 
holds. 

Remark 2. We recall that By [DS2], [DS3] if SM(£ ) holds for the potential (z,U) we can find a 
neighborhood O e {z) of z such that MUSM(O e (z)) holds for the interaction U. 

We stress that the above Remark 2 gives only a local implication. On the global side the 
relationship between MUSM(„4) and USM(„4) is not trivial. It is in fact possible to have a sort of 
layering phase transition which prevents MUSM(„4) to hold even though USM(„4) does hold. On 
the positive side, following an argument in the same spirit as [MOS] (i.e. that no phase transition 
may happen in a one-dimensional boundary of a regular two-dimensional domain), we rule out 
such a possibility in d = 2. We have in fact the following Theorem. 

Theorem 2.1. Letd = 2. Then USM(A) => MUSM(A). 

On the negative side we show that the aforementioned pathology may indeed happen. In 
Appendix A. 2 we give in fact an ad hoc example of an interaction U (in d = 3) such that: 

- GUSM holds; 

- there exist z and A of the form (2.10) such that the multi-grancanonical measure associated 
to (z, U) exhibits long range order. In particular there exist r, r' such that 



We finally mention that, in the context of the two-dimensional Ising model, the above Theo- 
rem implies the following. Consider a standard Ising model with a non-homogeneous external field 
which is however constant in cubes of side £; then for each (3 < (3 C , there exist I and L 3> £ such 
that SM(L) holds uniformly in the external field. 

2.8. Block Averaging Transformation (BAT) 

Let jj, z be the (unique) infinite volume Gibbs measure of a finite range translationally invariant 
lattice gas satisfying Condition SM(^o) and I an integer. In this case we can define the block 
averaging transformation directly in infinite volume. We partition the lattice L = 7L d into cubes of 
side £, i.e. C = [j i€Ce Qt{i). We recall that the random variable Ni has been defined in (2.12); it 
takes values in the set 




to? := {0,1,..., £ d } 



(2.14) 
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We then define the centered and renormalized random variables M, as in (1.3); it takes values in 
the set defined in (1.4). We finally let M_ := {Mi, i G Ce}- The BAT renormalized measure, 
that we denote by fJ,P , is then the (joint) probability distribution of M_ under fi z . Denoting the 
renormalized configuration by m = {rrii G &p ,i G Ci}, the measure fiP is formally given by 
IJ^z\ll\) = uAM = rn). We avoid the troublesome issue of describing Gibbs measures on non 

compact single spin space (see [EFS] for a discussion) and consider [ip only for finite £. Therefore 

- (£) 

the setup previously described applies to the finite single spin space Q] . 

It is also possible to use a finite volume setup. Given the integer p we will denote by A p CC C 
a cube with side 2d£p. We have A p = (Jig/ Qlif) where I p CC Ci is a cube of side 2dp. Let z be 
the finite volume Gibbs measure for our lattice gas with activity z enclosed in A p with r boundary 
condition. We denote by li) ' the finite volume renormalized measure arising from the application 
to // Z)P of the Block Averaging Transformation on scale £; it is defined as: 

*4''J dm, i e I p }) := MX P! , ({Mi = m z ,i£ I p }) , m, G ^ (2.15) 



2.9. Main results. 

We first discuss the case when the global Condition GMUSM holds. The most relevant example is 
the standard two-dimensional Ising model for T > T c . In such a case we are able to prove that the 
renormalized measure llP is, for each (finite) i large enough, Gibbsian w.r.t. a potential of 
bounded || • || Q norm (for suitable a > 0). We can furthermore control the i dependence of the norm 
II^^IU- We note (see [IS], [N]) that ^f 1 converges weakly to ^i^Ce^Pi where (fi denotes a standard 
Gaussian measure. Accordingly, should converge to the interaction of independent harmonic 
oscillators. Unfortunately, as the limiting interaction has not finite norm (since the limiting single 
spin space is unbounded), this convergence cannot be described in the || • || Q norm. However this 
lack of convergence affects only the (somehow trivial) short range part of the interaction; we will 
decompose the potential into a short and a long range part $^- ) = <I>W' sr + (3k G N : 

^W,sr _ q jjjgjj^jQ > K y ^-jjgjj introduce a large field cutoff (diverging as i — > oo) to control 
the short range part: it will converge to the potential of independent harmonic oscillator for values 
of the image variables within the cutoff. We note that this result would be false for large image 
variables. The precise statement is given in the following Theorem. 

Theorem 2.2. Let U satisfy GMUSM. Then there exists a > such that for any z G (0, oo) and 
t large enough fiP is a translationally invariant Gibbs measure w.r.t. a potential 3>^) for which 



< K{i) 



for some constant K(£) < oo. 

Furthermore it is possible to decompose the potential into a short and a long range part, 
$(t),sr + $(i),ir } such that 3k G N : <I>^ )!Sr = if diam(7) > k and the following holds: 
(i) for the same a as before 

lim <l' :,: - ,r - 
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(ii) there exists a constant a > such that 



lim sup ^P' sr (mj) 

\m I \<l a 



o, 



lim sup 



mi i<V 



(£),sr/ , 
{/} + 



for any / CC |/| > 2 



for any i £ £t 



When we assume only the local Condition MUSM(„4) our results are much weaker. Before 
discussing them, let us first note that for the standard two-dimensional Ising model this Condition 
holds for T <T C away from the phase coexistence line z = z* (z* corresponds to zero magnetic field 
in the spin variables), i.e. for each T < T c , MUSM(„4) holds for A = [0, oo) \ {z*}. We are not able 
to deal directly with the BAT defined in infinite volume, but we have to start from the finite volume 
transformation and take the thermodynamic limit. Moreover, we need to introduce a large field 
cutoff also in the long range part of the interaction. These difficulties are of course related to the 
limiting single spin space for which the usual (i.e. uniform in all possible b.c.) Gibbsian formalism 
do not apply. We refer to [EFS] for a discussion on the problems connected with the introduction 
of a norm for interactions defined on a non compact spaces. Our results are formulated as follows. 

Theorem 2.3. Let U satisfies Condition MUSM(A) and z > 0, z G A. Let also $( £ ' r ) be the 
(finite volume) potential associated to the (finite volume) renormalized measure l^i'j} ■ Then it is 
possible to decompose the potential into a short and a long range part, $( £ > T ) = <J>(^ T )> sr + <£>(^ T )> ir ; 



such that 3k G N : $ 



(i,T),sr 



if diam(I) > k and the following holds. There is a constant 



e = e(z) > such that for any I CC £i and any i large enough 



3 lim <l' : /- 7 ? ' •;//// • =: <!>/ (/„/) 



p— >oo 



uniformly for mj G Q) , |mj| < ey/ x\Qe\i T £ ^ 



3 lim $>\ l ' T ^' sr (mi) =: &\ l> ' b ' (mi) , uniformly for mj G £1^' , r G 



>(<) 



(2.16) 



Furthermore, there are a = a(z) > 0, a = a(z) > such that for the same e = e(z) as before 
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limVe^l sup <S>f' lr (mi) 

\mi\<ey/x\Qj 

<5> ( p' sr (mi) 



lim sup 

\mj\<e a 



lim sup 

* i. 



= 



\ mi \<i°- 



for any / CC £ e , \I\ > 2 



for any i G £,£ 



Warnings. 
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- Taking advantage of the symmetry of our Conditions w.r.t. the map z i-> we shall 
assume, without loss of generality, that all the activities are bounded by 1. This will be used 
extensively without further mention. 

- We denote by C a generic positive constant whose numerical value can change from line to 
line. From the statements it will appear clear from which parameters it depends on. 

3. Computing the renormalized potential via cluster expansion 

In this section we discuss the BAT transformation in finite volume. We will compute the renor- 
malized interaction via a cluster expansion: the convergence of the expansion will be ensured by 
the validity of condition Cl(m, £(d)) for the constrained (multi-canonical) systems. This condition 
CI, in turn, will be deduced from the MUSM property of the original system in Section 5. 

To simplify notation we write the Boltzmann's factor (with r boundary condition) for a 
configuration rj in the volume A, 77 G {0, 1} A , as exp (+Ha(t]\t)) where 

H A ( V \r) := -E A (r,\r) (3.1) 

Let us set L := d£; given the odd integer p, let A p be the cube with side 2d£p given by 

' { x = (xi , . . . , x d ) G C : -d£ (p + |) + dl + 1 < Xj < +di (p + |) , j = 1, . . . , d} d£ even 

k {x = (x lt ...,x d ) €£:- {d£p+^-)+d£ < Xj < d£p+^-,j = l,...,d} d£ odd 

We can write A p = U i€ i p Qe(i) where / p CC £{ is the cube of side 2dp given by 

' {i G Ci : —d (p + i) + d + 1 < Xj < +d {p + |) , j = 1, . . . , d} if d is even 
k {i G C t : - {dp + + d < xj < dp + 4=1, j = 1, . . . , d} if d is odd 
Let us introduce the quantity: 

where n = {n i; i G I p } G Q { P := O ie/p fif ) = {0, 1, . . . , £ d } Ip and 



A p 



I P = < 



ni) :={r?G{0,l}^W : N^r,) = m] , fig> := (g) SI™ (3.3) 



It is convenient to look at the renormalized measure fJ^ T J in (2.15) in terms of the variables n; 
such measure is Gibbs w.r.t. to the renormalized Hamiltonian given by 

Hfj\n) = log Zfjl (3.4) 

Given n G fi^; we can look at the quantity Z^'J^ defined in (3.2) as the partition function of 
a (generally not translationally invariant) system which is the original lattice gas constrained to 
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have fixed values of the total number of particles in each block Qe(i), i G I p . Its elementary 
configurational variables are the original spin configurations in each block Qe(i) compatible with 
the assigned value rii of N{ namely the set njp^ defined in (3.3). The elements of f^™ 1 - 1 will be called 
block variables not to be confused with the renormalized variables rtj. We also call multi-canonical 
these constrained systems. 

We will consider blocks of these block variables of size d; these corresponds to the blocks Qh{i) 
with L = di in the original variables. The reason for this choice will appear clear in the following 
sections: it corresponds to the minimal size for which we are able to prove, for the constrained 
model, the validity of our Condition Cl(m,e(d)). In other words, to meet Condition Cl(m,e(d)) 
we have to choose m = d and I sufficiently large. With respect to the general setting of Section 2 
we have Jlj = £l\ whereas the potential is the one inherited by the original model. In particular, 
if we choose i larger that the range r of the original interaction, then only contiguous blocks will 
interact. We repeat that the size of the blocks that in Section 2 was generically called m now equals 
d. The main result of this section is stated as follows, where, for V G P^ k \i), we let V CC £g be 
such that V = U jev Qe(j)- 

Theorem 3.1. Consider a d-dimensional lattice gas with finite range, translationally invariant 
interaction. Let I G N and suppose there exists a closed V C [0, 1] such that 



sup sup sup sup sup 

ielp k=i,-,d VeP w^ ReT) w <r,C,r 



Z v , R (a^\a^-),r)Z VtR (C^\C^-\r) 
Z v ,n {a(^+),C (k '-Kr) Z V ^ (C( fc >+\ a(*.->, r) 



(3.5) 



where vf := (\Qi\V) v n and 5{i) -> as t -> oo. 

Then, the measure [if' r) z defined in (2.15) is Gibbsian w.r.t. a potential $C> T ) = {$^' r) , X C I p }. 
Let 



x 



x 



We have the following: 

(i) For each X C I p with dt(X, 1^) > 3d and mx G ' , ^x'^ does not depend on r (and I p ). 
In particular for each Icc£< 



3 lim ®x' T \m x ) =■ ^x'( m x), uniformly for m x G M x ' -,t G f2 



V) 



p— >oo 



(3.6) 



(%%) Let&W = X CC C e }, we have a decomposition = $(<)>«•+$(<).''■ where $^ ),sr = 

if diam^(X) > 3d and there are constants a > 0, C such that: 



E 



e a l x l sup 



X30 m x €M 



<$> { x ] ' lr {mx) <CS(£) 



(3.7) 



Remark. The potential <J>^' ^ will be explicitly constructed (see (3.70) and (3.69)) below. In Section 
5 we will show that we can take $^' sr (mj) = —ml/2 and there exists a constant a > such that 
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for each X, \X\ > 2 



sup 

rn x eM { £> 
|m x |<*" 



with j(£) -^Oifi^oo. 



Similarly to what has been done in [HK], in order to compute the renormalized potential 
and prove Theorem 3.1, we are going to apply to the constrained systems the method developed 
in [O], [OP]. To simplify the exposition we will treat in detail only the two-dimensional case. An 
analogous treatment can be developed for the d-dimensional case along the lines of [OP]. For the 
same reason, we discuss only the case of periodic boundary condition in A p ; the case of general 
boundary condition can be treated along the same lines with minor changes giving rise to estimates 
uniform in r. 

In the rest of this section we will express the coordinates of points and components of vectors 
in Cl in L units. Let us denote by ei, e 2 , respectively, the coordinate unit vectors in Cl'- &i = (1,0) 
horizontal, e 2 = (0, 1) vertical. Recall that since now d = 2, we have L = 21. We further partition 
Cl into four sub-lattices of spacing 2L: 

C-L = C-2L U £$L U £%L U ^2L 



where: 



L '2L 
L '2L 

c c 

Z -2L 



L -2L 



{i = (i 1 ,i 2 ) G C L :i 1 = 2j 1 ,i 2 = 2j 2 ,for some integers ji,j 2 } 

(3.8) 



= C4l + e 2 



C£ L + ei + e 2 = C% L + ei 

+ ei = £% L + e 2 = Cf L + ei + e 2 

We also set, for i G £l: 



i4 i :=Q L (2i), 5, :=Q L (2i + e 2 ) C, := Ql(2i + ei + e 2 ) A := Ql(2« + ei). (3.9) 
(See Fig. 1). 

Then we can partition the torus A p into the union of the L-blocks of the four types: A, B, C, D: 

A p = A p U B p U C p U V p 

where 

A p :={Ai : \ij\ < J = 1,2} 

and similarly for Bp,C p ,V p . 

Given a renormalized configuration of our multi-canonical model and a block Ai we denote by 
ai a generic original lattice gas configuration compatible with the four renormalized configurations 
{rij, j G £e : QKj) C ^4^}; in other words: a, G ( X>j : Q f (j)cA i ^j n3 ' > (recall that £ = -|). Similarly 
for /3j,7j,<5j. We simply denote by a, (3,^,5 the configurations in Ap, C p , P p , respectively. 
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Fig. 1 

Let us now quickly describe our strategy. We want to transform the constrained system into 
a polymer system (see, for instance, [GrK], [KP], [D3]) which, by condition (3.5) will turn out to 
be in the small activity region. To be more precise we shall prove the following formula: 

4\ = 4\4\ (3.io) 

~ (£) 

where n is factorized in the sense that it has the form of a product of partition functions 
in suitable volumes not depending on p; the dependence on n of the single factors is local. The 
partition function Z K ^ n describes the reference system around which we perform a perturbative 

(£) 

expansion. On the other hand, Er A J n is the partition function of a gas of polymers; it has the form 



C = l + E E UCnM) (3-11) 

fc>l Ri,—,R k j=l 



where the polymers Rj, that will be defined below, are geometrical local objects living on scale 
L = 1t\ the sum in (3.11) is restricted to "non-intersecting" polymers so that the unique interaction 
between polymers is a pairwise hard core exclusion. Finally the activity Cr 3 (?i) depends only on 
the rtj's with i localized on the polymer. It is already clear from this preliminary discussion that 
expression (3.11) is well suited to compute renormalized potential: in order to get good estimates of 
the norm of the renormalized potential we shall need that the polymer system described by n 
is in the small activity region. 

To get expression (3.10) we will perform a sequence of block decimations like in [O], [OP]. 
We start by integrating over the 5-variables, then the 7-variables, the /3-variables and, finally, the 
a-variables. Using Condition (3.5) we will be able to prove that at each step of decimation the 
system described by the surviving variables is weakly coupled. 
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We use the following notation for the interaction (which is defined independently of the 
multi-canonical constraints) between two sets Ai and A 2 : 

W KlM (r} Kl \r} K2 ) :=W(r] Al \r] A2 ) = H AlUA2 {va 1 , Va 2 ) - H Al {rj Al ) - H A2 (r] A2 ) (3.12) 

where rj Al ,r) A2 € {0, l} Al , {0, 1} A2 , respectively. Recalling that L is larger than the range of the 
interaction, we can write: 

H Ap (a Ap )= H A kl M+ E H Bk2 {(3 k2 ) + W Bk2>Ap {(5 k2 \a) 

ki:A ki eA p k 2 :B k2 €t3 p 

+ E H c k3 (7k 3 ) + W Ck3 ,A P uB p {lk 3 \f3,a) ^ 313 ^ 

ks:Ck 3 dC p 

+ E H D H (.S k4 ) + W Dki!ApUBpU c p (S k Al,P,a) 

fc 4 :Dfc 4 eE> p 

Again the above decomposition of H holds independently of the constraints on the number of 
particles in the blocks; in (3.13) we have only used that L > r so that there is no direct interaction 
between blocks belonging to the same sub-lattice. 

To simplify notation we shall often omit from H, W the subscripts referring to the various 
domains; the symbols used for the arguments of the functions H, W should be sufficiently clarifying; 
moreover we will also omit the explicit extensions of the sums (or products) over k\ , k 2 , k 3 , fc 4 as 
well as the one over a 6 <S>i : Q f (i)cA p ^ ni \ an d similarly for (3,J,5. We have: 

Z K,n = E II ex P )) E II ex P ( H ^ 2 ) + W((3 k2 \a)) 

Q ki P k 2 

(3.14) 

x E II ex P fas ) + I P> <*)) E II ex P ( H ) + W(5 k4 1 7 , P, a)) 

7 fc 3 <5 fc 4 

We first perform the sum over 6 variables; using that the sums over different 5 ki are decoupled 
since the size L of the blocks is larger than the range of the interaction, we get: 

< E E E II ((& 7) u . w, i)\ «)) (3.i5) 

a. p 7 ki 

where by Zjj k ((/?, 7)", (/3, 7)^, a)) we denote the partition function in D ki with boundary con- 
ditions (/3, 7)" on the top (up) and (/3, 7) d on the bottom (down) of D ki (see Fig. 2). More 
explicitly (/3,j) u is given by the restriction of /?, 7 to (simply called configuration in): Qh{2k^ + 
e2)U(5L(2fc4 + e2 + ei)UQL(2fe 4 + e2 + 2ei) = B k4 UC k4 UB k4+ei whereas (/3,^) d is the configuration 
in Q L {2k A - e 2 )\jQ L {2k A - e 2 + ei) U Q L {2k A - e 2 + 2e l ) = B ki _ e2 UC fc4 _ e2 UB ki+ei _ e2 . Finally 
a in (3.15) denotes the configuration in (5l(2/c 4 ) U Ql(2/c 4 + 2ei) = Afc 4 U A k4+ei . 

Notice that we are also presently omitting the explicit dependence on n and L. Let denote 
a given reference configuration in fi^. We write 

V 7 (3.16) 

Z Dki {(^r,(0y,a))Z Dki ((0r,(/3, 7 ) rf ,a)) 

Z D . ((())«, (0)««, a)) 
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Where by (0) u ,(0) d ,a we mean the boundary condition on D ki obtained from ((3,j) u ,((3,j) d ,a 
by substituting (f3, 7) with (0) both in the "up" and "down" blocks; similarly (/?, 7)", (0) d , a ; 
(0) u , (/?, 7) ,a denote the boundary conditions on D ki obtained from {/3, r y) u , (f3,j) d ,a by substi- 
tuting ([3, 7) with (0) only in the "down" , "up" blocks, respectively. We call the above operation 
"splitting" of the partition function Zo ki ((/?, 7)", (f3, ^) d ,a)) in the vertical e2 direction. We set 

*<*> ( R*. = ZD ^ ((^r,(Pn) d ,a))Z Dki ((Or, (Of, a)) 
D( ' P,7j ' ^ fc4 ((/3,7)-,(0)^a))Z Dt4 ((0)-,(/3, 7 ) d ,a)) ^ ° 

The quantity (a, 7) can be considered as an effective interaction potential between a,/3, 7 
variables coming from decimation of the 5 variables. In what follows we will exploit condition (3.5) 
above to deduce that <E>^ (a, (3, 7) and other similar quantities are uniformly small. 
We can write: 

Z K = E • • • E • • • E II ex P Wto ) + W( lks a)) 

a 7 k 3 

x Z D , 3+e2 ((0r,(/3,7 fe3 ) d ,«) Z D , 3 ((/3, 7fc3 )«,(0) d ,a) (3.18) 
x [] [Z Dki ((0)-, (Of, a))] + < 4 («, A 7)) 

In (3.18) above we associated to every C k3 block in C p the two terms Zo k3+e2 ((0)", (j3, J k3 ) d , a) , 
■^-Dfc 3 ((/?, 7fc 3 ) U ) (0) d ,a) coming from the splitting of the original partition functions over the vol- 
umes Dk 3 +e 2 ,Dk s , respectively. Notice that 

J>xp (H( lk3 ) + W( lk3 %a)) Z Dks+e2 ((Or, ((3,lk 3 ) d ,,a) Z Dks (((3,^ k3 ) u , (0) d ,«) 

7fc s (3.19) 
= Z 5k3 ((0) u ,(0) d ,a,(3)^Z s ^ ((Or,(O) d ,(a,(3) l ,(a,0r) 

Where: 

C k3 := C k3 U D k3 U L> fc3+e2 (3.20) 
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is a 3L x L rectangle DCD centered at C k3 (see Fig 3) and Z~ ((0) u , (0) d , (a, f3) l ,(a, f3) r ) is the 



partition function in Ck 3 with (0) boundary condition on the top and on the bottom; (a, j3) 1 on the 
left and (a, (5) r on the right of Ck 3 - Here by "on the left" of Ck 3 we mean "in A k3+e2 U B^ 3 U Ak 3 " 
and by "on the right of Ck 3 we mean "in Ak 3 + e2 +e 1 U-Bfc 3+ei UA)- 3+ei "; see Fig. 3. In what follows 
we will continue to use "on the top", "on the bottom", "on the left" and "on the right" for the 
boundary conditions to a volume in a similar sense. 
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Dk 3 
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Fig. 3 

The operation described by equation (3.19) above is called "gluing" of the partition functions 
Z D H+e2 {(0) u ,((3,lk 3 ) d ,a), Z Dk3 (([3,~f k3 ) u ,(0) d ,a) on C k3 in the vertical e 2 direction. 
Now if in (3.18) we multiply and divide by 

[]% fc3 ((0r,(0) d ,(a,/3) z ,(a,/?r), 

we get: 

z Zn = e n ex p ( h m) e n ex p ( H tf>* ) + w (^ 2 \a)) n z ^ k3 w> (°) d > («> ^ («> 

a fcl /3 fc2 fc 3 

e /v^t) n (! + < 4 («. & 7)) n («))] _1 

7 fc 4 fc4 

(3.21) 

where (7) is the product (Bernoulli) probability measure on 7 parametrically depending on 
a, (3 given by: 

/^(7):=n^f 3 (7 fc3 ) (3-22) 

fc 3 
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where 

< fe - ): = Z5ij ((o)., ( o)<«!( a , /3 y,( a , f3 ) r ) exp(ff( ^ ) + w ' < *> l ' 3 - a)) (323) 

At this moment we operate again a "splitting" but now we act on the partition function Z~ ((0) u , (0) d , (a, (3) 1 , (a, f: 

3 

in the horizontal e\ direction; namely we write: 

((0) u ,(0) d ,(a,f3) l ,(a,f3Y) 

'Z^ {(Or,(0) d ,(a,py,(a,PY)Z 5k3 ((0)" , (0) d , (0) 1 , (0) r ) 

z~J 3 ((o)^(o)^(a,/3)^(o)-)z~ 3 (( 3 o)^(o)^(o)^(a,/3r) " L + 1 ; cud 

Z^ {(0r,(0) d ,(a,P) l ,(0y)Z 5k3 {(0r,(0) d ,(0) l ,(a,py) 



X 



z 5k3 ((oy,(o) d ,(oy,(oy) 



We set 



Z~ {(Oy,(0) d ,(a,(3) l ,(a,(3y)Z~ ((Of, (0) d , (0)< , (Of) . 

z 5k3 ((oy,(oy,(a,py,(oy)z 5k3 ((oy,(o) d ,(oy,( a ,py) ' ' 1 ~ ' 

We remark that: 



^2exp(H((3 k2 ) + W((3 k2 \a))xZ 5k2 ^ ((0) u ,(0) d ,(0) l ,(a,(3 k2 ) r ) 
% ((0)-,(0) d ,(a,^ a )',(0r)=Z~ ((0),a) 



(3.26) 



where B k2 is the set, centered at B ka , having the shape of a capital H given by: 

Bk 2 '■= B k2 U Cfc 2 U C k2 - ei U D k2 U D k2 _ ei U Z)fc 2 _ ei+e2 U D k2+e2 (3.27) 

see Fig 4. The above operation, described in (3.26) above, is a "gluing" of the partition functions 
Z~ ((0r,(0) d ,(0y,(a,p k2 ) r ), Z~ ((0y,(0) d ,(a,(3 k2 y,(0y) onB k2 inth eei direction. 

The boundary condition on B k2 in the partition function Z~ ((0), a) are everywhere except 
for the A-blocks A/ C2+e2 , A k2 touching on the top and on the bottom, respectively, the block B ks . 
We write: 

Z~J{Q),a)=-.Z~ k ^),{a) u M) d ) (3-28) 

with (a) u , (a) d , given, respectively, by the restriction of a to A k2+e2 , A k2 . 

Similarly for the term Z Dk4 ((0)", (0) d , a)) appearing (at the power —1) in (3.18) we can write 

Z Dk4 ((Or, (0) d , a)) = Z Dk4 ((0y, (0) d , (a) 1 , (ay)) = Z Dk4 ((0), (a) 1 , (a) r )) (3.29) 
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Fig. 4 

where by (0), (a) 1 , (a) r we mean the boundary conditions, outside D^ A given by everywhere 
except for the two blocks Ak 4 , A/- i+ei , contiguous to Dk 4 ; (a) l ,(a) r are the restrictions of a to 
A ki , A ki+ei , respectively. 

Now we perform a "splitting" in the e\ direction of the quantity [Zd^ ((0), (a) 1 , (a) r ))] ; 
namely we write: 



[7 (m( y r . n x r i ( zd* ((o),«(on)^ 4 ((o),(oy,(«r)) 

[Z Dfc4 ((0),( a) ,(a)))] =[ ZDti(( o) i(a)I)(an)K((0)i( o ) . (on )- 1 + 1 y 

z DkA (( ),(oy,(on) 



(3.30) 



Z Dki ((0),(ay,(0r))Z Dt4 ((0),(0)', (a)')) 

We set: 



Z Dfc4 ((())»', (OP) Z g>4 ((0),(0)',(ar)) \ (4) 
Z D . ((0), (a)', (a)')) Z Dfc f(0),(0)<, (O)-)) 



We introduce the product probability measure [1% ((3) on (3, parametrically dependent on a, 
given by: 

^(P)--=H» a B k2 (Pk 2 ) (3.32) 



where 



»K &>») := Z~ ((0),(a)",(a)') 6XP + 

xZ~ 2 _ ei ((0r,(0) d ,(0)'(a,/3 fc2 r,)Z~ 2 ((0r,(0) d ,(a,/3 fc2 )',(0r) 
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Now we proceed similarly to the step leading to (3.21). We multiply and divide the expression 
on the r.h.s. of (3.21) by 

\{Z~ k ((0), («)«,«); 



k 2 



by inserting in the r.h.s. of (3.21) the expression given by (3.33) and after operating the splitting 
described in (3.24), the gluing described in (3.26) and the splitting described in (3.30), we get: 



z A p ) ,n=EII ex p(^ 1 )) m,M l ,(oy))z Dki ((o),(oy,(a fcl r))] _1 

a ki 

x u z B k2 ((o), x>?(/*) n [ z c H ((o))] _1 n i 1 + *§,) e^'^t) 

k 2 /3 k 3 k 3 7 

n (i + < («, a 7)) n + < («)) n ^ «°))] _1 



X 

(3.34) 

where we used the shorthand notation Z~^ ((0)) for Z~^ ((0) u , (0) d , (0)', (0) r ) and Z Dki ((0))) for 

z Dki ((or,(o) d ,(oy,(oy). 

Now we perform, on the partition function Z~ ((0), (a) u , (a) d ), a splitting a bit different 
with respect to the previous ones. Let F^ 2 be the horizontal L x 3L rectangle CBC contained in 

Fk 2 = B k2 U Ck 2 U Cfc 2 _ ei (3.35) 

We can write: 

Z~ h 2 ((0), (a)«, (a) d ) = £ exp + (tf ((<5) d ) + W((a)«|(*T)) + ^((a) d |(^))] 

(<5)",(<5) d 

( Z Fk2 ((0),(a,5r,(a,5) d )Z Fk2 ((0), (0)" , (0) d ) \ 
X ^ Ffc2 ((0),(a,^,(0)^^ fc2 ((0),(0)«,(a,^) + ) 

Z Fk2 ((0),(a,Sr,(0) d )Z Fk2 ((0),(0)",(q,^) 
Z^((0), (())«, (0)««) 

(3.36) 

Where, for a generic 5 G ®j-.Q t {j)cv, we denote by (5)" the restriction of <5 to Dfc 2 _ ei+e2 U 

Dk 2 +e 2 whereas we denote by (5) d the restriction of (5 to Df, 2 - ei U .Dfc 2 ; by (0), (a, (5)", (a,(5) d we 
mean boundary conditions on given by (a, 5) u on the top, (a, 5) d on the bottom and elsewhere 
(see Fig. 5). 

Let , F^ be the "horseshoe" shaped domains given by: 



F k? '■— B k 2 U Ck 2 U Ck 2 - ei U Z?fc 2 - ei +e 2 U D k2 + e2 , 

(d) (3-37) 
F k2 :— i?fc 2 U Ck 2 U Cfc 2 _ ei U -Dfc 2 U D k2 _ ei 
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Fig. 5 

(see Fig (6)). From (3.36) we easily get: 

Z F( „,((0),(a)«)Z F(d) ((0),(a) d ) 

Z~ ((0),(a) u ,(a) d ) = k2 - w " 2 , , (1 + *b (a)) (3-38) 

s fc2 W /'V 7 v. ; ; Z Fk2 ((0),(0) u ,(0) d ) V B *a v V v ; 

where Z p ( U ) ((0), (a) u ) is the partition function on the domain with boundary conditions 

everywhere except for Ak 2 + e2 where they take the value (a) u ( = the restriction of a to Ak 2 + e2 )', 
similarly Z p (d) ((0), (a) d ) is the partition function on the domain with boundary conditions 

everywhere except for A^ 2 where they take the value (a) d ( = the restriction of a to Ak 2 ); finally 
®bI ( a ) i s defined as: 

= < 5 £ s)a ,W ' (^((0),(a,fl-.(0WZ ft ,((0),(0).,(a.W " 'J 

where Jif ((5) u ^(5) d ) is a probability measure on <8).. /-^t> or parametrically dependent 

on a fc2+e2 ,a fe2 given by: 

W) =^P + (^(W") + W((ar\ (5)")) + W((«) d |(5) d ))) 

Z Ft2 ((0),(q^)",(0) d )Z Ffc2 ((0),(0)",(M) d ) (3.40) 
X Z F( „)((0)>)")Z F(d) ((0),(an 

Indeed /Ik 2 ((<5) u , (8) d ) has the form of a product measure over the "up" and "down" variables but 
in (3.39) we are averaging, with respect to , a function which couples these variables so that the 



result is a & Bk (a) which is a non-factorized function of (a) u , (a) d . 



s \ IY ) W 1 Ul II i (I J..IA '.I..I UH I v n 1 /.. V L I I 
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Fig. 6 

By inserting (3.38) into (3.34) we get 

z A p ) ,n=EII ex P(^ 1 )) ((o),M l ,W))z Dki . ei ((oMo^K.n)]" 1 
z fM ((o),(« tl )")V' (W^KJ^n^ + ^W) [^ 2 (o)] _1 

E ^(/?) II [% fcs ((0))] _1 II (l + < } 3 (a, /?)) 

P k 3 k 3 

e n i 1 + < 4 («. ^ ^)) n ( x + *z w) n ^ «°)) 



where we have used the shorthand forms ((0), (a^)')), respectively Zr> k _ ei ((0), (afcj r )), for 
Z Dfci ((0),(a fcl )',(0r)), ^D fcl - ei ((0),(0) I ,(a fci r)) and Z Fk2 (0) for Z Ffc2 ((0), (0)«, (0) d ) . 

We notice that if in (3.41) above we neglect all the "small quantities" <3? and *S> and we 
use that n^iP) an d M3 (7) are normalized measures, then, by performing the sum over the 7,/? 
variables, we get a factorized partition function describing a system of independent a variables. 
So we substantially have already reached our goal; we want now to manipulate a little bit these 
factorized terms (the product over ki) in order to get a simpler expression with a more transparent 
physical meaning. 
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We use the notation A kl to denote the 3L x 3L cube centered at the block A kl : 

A kl := Q 3 l (2fci), fci G £ L (3.42) 
Let G kl denote the annulus obtained from by removing the block itself: 

Gfcj := A kl \ A kl = B kl U B kl _ e2 U Gfcj U C kl - ei U G^-^-^ U Cfc 1 _ e2 U D kl U Ufcj-ei (3.43) 

We denote by Zc; fcl ((0),afe 1 ) the partition function on G^ with boundary conditions a kl 
on the "hole" A kl and elsewhere. Moreover let Zr> k u£> k ((0), a^i , {fil) w ■, {fll) d ) denote the 
partition function on the (non-connected) set D kl _ ei U D kl with boundary conditions a kl on A kl , 
(/3 7 ) u on the up part of G kl \ (D kl _ ei U D kl ) (namely in G fcl _ ei U B kl U G fel ), (/? 7 ) d in the down 
part C kl - ei - e2 U -Bfci-e 2 U Ck 1 - e2 and elsewhere, (see Fig (7)). 



(A 7) 



(ft 7) 



Fig. 7 

Indeed we have the following factorization: 

Z Dkl _ ei uD kl ((0), a kl , (/3 7 r, (/3 7 ) d ) = Z Dk± _ ei ((0), a fcl , (/3 7 ) u , (/? 7 ) d ) ((0), a fcl , (/? 7 ) u , (/3 7 ) d ) 

(3.44) 

We have: 

Z Gfcl ((0),a fcl )= ^ exp(i7((/? 7 D+i7((/? 7 ) d ))Z Dfci _ eiUDfci ((0), a fcl , (/3 7 ) M , (/3 7 ) d ) 

(/3 7 )",(/3 7 )^ 

(3.45) 
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We can write: 



z Gkl ((o), akl ) = Yl ex p ((/w + mhf)) 

(/3 7 )",(/3 7 )rf 

^D fcl - ei uD fcl ((0), » kl , (£7)", (/?7) d ) Z Dki _ eiUDki ((0), « fcl , (0)", (0) d ] 
Z Dkl _ ei uD kl ((0), a fcl , {h) u , W) Z Dki _ ei yj Dki ((0), a fcl , (0)«, (/3 7 ) d ) 

((0), a fcl , (/3 7 ) M , (0) d ) ^ Dtl _ ei uD tl ((0), a fcl , (0)«, /3 7 ) d ) 
((0),a fcl ,(o)",(0)<*) 

Z F( „)((0),(a fci r)Z F(<i) ((0),(a fcl ) d ) 

^ ^2 / 1 , &{]) ( \ ) 

z Dki _ eiUDki m,* kl ,(or,m v + ^ a *>) 



where 



(/3 7 )-,(/3 7 ) d 

1 uD fcl ((0),a fcl , (/3 7 ) u , (/3 7 ) d ) 2 Dti _ eiUDfei ((0),a fcl , (0)", (0) d ) 
Z Dkl _ ei un kl ((0), a fcl , (/J 7 )«, (0) d ) Z Dfci _ eiUDfci ((0), a fcl , (0)«, (/3 7 ) d ) 



and 



Ml? 1 j ((/37)", (/?7) d ) := exp (#((/? 7 D + # ((/3 7 ) d )) 

x ^ 1 _. 1 up tl ((0),« fcl ,(/3 7 )",(0) d )Z Dfci _ eiUDfci ((0),a fcl ,(0)",(/37) d ) 
Z F( „)((0),(a fcl )-)Z F(d) ((0),(a fcl ) rf ) 



We write: 



From (3.46), (3.47), (3.48), (3.49) we get 



Z F( „)((0),(a fci r)Z F(£i) ((0),(a fcl ) d ) 



We define the Bernoulli probability measure n\ (a) as 

Mi (") := n^i(° fc i) 



where 



^((0)) 



exp(H Ahi (a kx )Z Ghi ((0),a kl ) 
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in which by Z~ ((0)) we denote the partition function in with boundary conditions. 
In conclusion, from (3.41), (3.42), (3.50), (3.51) and (3.52) we get : 

z aL = n z A kl ((o)) n (o)] _1 n [% fc3 «o»] _1 n ^ «o» 

" 1 &2 ^3 &4 

e a*i («) n (i + )) n (i + < («)) n + *s 4 w) (3.53) 



X 



fci fc 2 fc4 



e ^ 2 ( a 2 ) n + < } 3 (a, a) e ^(7) n (! + < 4 («. & 7)) 



/3 fc 3 7 ^4 



We write 



with 



^Ap,™ — 4p,n"Ap,n (3.54) 

z aL ■■= n ((o)) n ^ _1 n [% fc3 «°))] _1 n ^ cw) (3.55) 

fci fc 2 k 3 fc 4 

and 

*2* -E»wn('+ < <■*. >) n (» + < <<->) n (» + < («>) 

a fci fc 2 fc4 

e ^ 2 (aj n { i + *§, («, /?)) e ^'"(7) n + < 4 («. a -r)) 



(3.56) 



x 

' P fc 3 



We are now ready to express n as the partition function of a gas of polymers whose only 
interaction is a hard core exclusion. 

We have to analyze the various interaction terms (the <£'s and ^'s ) appearing in (3.56). 
We see from (3.17) that the term (a,/3, 7), involving the a, (3, 7 variables in the annulus 
Q3l(2A;4 + ei) \Dfc 4 , corresponds to an "eight body" interaction among the A, B, C blocks adjacent 
to Dk A \ we see from (3.25) that («,/?) is a six body interaction involving the A and B blocks 

adjacent to Cfc 3 ; (a), (a) are two body terms involving the pair of A blocks contiguous 

to Bk 2 , -Dfc 4 , respectively. Finally («fci) is just a one body term. 

ki 

Looking at (3.22), (3.23) we can say that <J>^ (a, (3, 7) extends its action to all A and B blocks 
adjacent to the C blocks in Q^L^k^ + e\) (see Fig. 8), becoming a "twelve body" interaction. 
Indeed we have to average <3?^ (a, f3, 7) with respect to the product of the measures /x^ 3 (7fc 3 ) 
which are parametrically dependent on the a, (3 variables adjacent to Ck 3 ■ On the other hand, 
looking at (3.33), it is easily seen that we do not have to extend any more the region of influence 
of $^ (a, 7) because of the parametric dependence on a of (flk 2 )- Moreover, still looking 



at (3.33), we easily see that also the term (a, (3) does not extend at all its influence. Of course 



Ofc 3 

<3?£k (a), ^Dfc ( a ) ^aI ( a ki) do not extend, as well, their action. 
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Fig. 8 



So it is natural to define different kind of (many body) bonds corresponding to the above 
interaction terms. As a consequence of the above discussion we have the following kind of bonds; the 
bond , to which corresponds the weight <E>^ (a, (3, 7) which is given by the set of A, B and C 
blocks contiguous to Dk 4 united with the other A blocks adjacent from the exterior of (J3l(2/c4 + e\) 
to the already considered B blocks. So a Dj^-bond contains twelve blocks. We similarly define 

(now without any extension) the bond with weight ((3, a); the bond B^ with weight 
^bI ( a )' t ne bond with weight (a) and the bond A^ with weight ( a kt)- 

Given a bond b of one of the above kinds we define its support b as the subset of C obtained 
as the union of the Ql blocks making part of b. For any bond b we denote by the corresponding 
weight. Notice that will be, in general, a function of the a, (3, 7 variables associated to the blocks 
in b. For instance a bond b = can be seen as an element of (£l) 12 whereas b is a subset of 

the original lattice C given by the union of the twelve interacting blocks. 

We say that two bonds 61,62 are connected if b\ n 6 2 / 0. A polymer R is a set of bonds 
6fc which is connected in the sense that V i, j: 1 < i < j < k there exists a chain of 
connected bonds in R joining 6j to bj namely 3 6^, . . . bi h ; b im 6 R, m = 1, . . . h, 6^ = 6j, 6 ih = 6^: 
6- n bi / 0, m = 1, . . ./t- 1. 

The support i? of a polymers i? = 61, . . . , is simply i? = U^ =1 6j. We call 7£a p the set of all 
possible polymer with support in A p and 1Z the set of all possible polymers with arbitrary support 
in C. Two polymers i?i,i?2 are said to be compatible if R\ n R2 = 0; otherwise they are called 
incompatible. 
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Given a polymer 72 = 61, . . . , b& we define its activity £r as: 



k 

Cr ■= E^^E^^E^'^^II^KAT) (3.57) 

a [3 7 i=l 

Notice that, due to the Bernoulli character of the above probability measures, we can, as well, 
write: 

E M lfS («s) E E ^(Tr) I! & («s» ^5. -rs) ( 3 - 58 ) 

a~ B~ 7~ j=l 

R R R 

where a~, /3~,7~ denote the a, /?, 7 variables in i?; 

a~,/3- 

K (7) = 



\i*W= II "c£(tO (3-59) 



and so on. 

Going back to the specific structure of our multi-canonical model it is immediately seen that 
the activity of a polymer R is a function of the renormalized variables n, (= number of particles 
fixing the constraint in the block Ql(i)) only for Qh{i) £ R- To make explicit this dependence we 
write 

Cr = ( R (.n~) (3.60) 

where n~ = {ni} 0t(i ) C £- 

From (3.56), (3.58) we get the desired expression: 

k 

7W 



-R i nR 3 -=0,i< 3 = l,...,n 

Now we state a Proposition referring to a general class of polymer systems. Its proof, which 
is based on the standard methods of the theory of the cluster expansion, can be found in [O] (see 
also [GMM], [KP], [D3], [NOZ]). 

Proposition 3.2. Consider a general polymer system (see [GrK], [KP], [D]) where the only in- 
teraction is a hard core exclusion forbidding overlap of the supports R of the polymers R. Its 
partition function is: 

k 

h a = 1 + E E^ lie* ( 3 - 62 ) 

fe ^ 1 ^ Hi J ...,H fc :H i cA, 1=1 
R i nRj=H,i<j = l,...,n 

Suppose that: 

i) 3 k > such that the number of different polymers R with m bonds (we write \R\ = m) and 
support R containing a fixed point (say the origin) is bounded by K m ; 
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ii) 3s >0 such that \( R \ < e^L 
Let 

<p T (Ri,---,Rn) = - } (-l) #edgesin9 (3.63) 



g eG{Ri,...,R n ) 



where G(R\, . . . , R n ) is the set of connected graphs with n vertices (1, . . . ,n) and edges i,j 
corresponding to pairs Ri,Rj such that Ri Pi Rj / (we set the sum equal to zero if G is 
empty and one if n = 1). If 



1 x 
e < — e 

K 1 + X 



x=(52-I) 



(3.64) 



then there exists a positive constant C(e) such that 



\R\ 



£ \MRi, ■ ■ ■ ,Rn)\ f[ IC*I < C(e) (e exp j ^L_i j, ) ( 3 .65) 



R 1 ,...,H„:R i CA 1 1 



-A = exp < 



£ MRi,---,Rn)U^\ ( 3 - 66 ) 



^ -Rl , ■ ■ - , Rn '-Ri CA 



In our context, it is clear that we can find a constant k so that the hypothesis (i) of Proposition 
3.2 holds. It is also clear from (3.17), (3.25), (3.31), (3.39), (3.47), (3.49) that there exists a universal 
constant C such that hypothesis (ii) holds with e = C5(£) (recall that 5(1) — ► as ^ — > oo) so 
that (3.64) holds for any ^ sufficiently large. In fact in the two-dimensional case we use a weaker 
condition: we do not need, in the left hand side of (3.5) to take the supremum over V £ P^\j), 
but only the analogous condition only for the squares Ql and for the L x 3L rectangles. 

Then, using the results of Proposition 3.2, we can compute the renormalized potential and 
perform the thermodynamic limit. Suppose, instead of considering periodic boundary conditions, 
we had a generic b.c. r outside our cube A p . It is clear that we can apply the same procedure 
(block decimation and cluster expansion) that we have used above in the case of periodic boundary 
conditions and get very similar results. Let us briefly sketch the differences. 

Recall that our square A p has a side being an integer multiple of the elementary square Ql 
with side L = 21; then certainly we will have a horizontal edge of dA p adjacent to a row (of thickness 
L) made by C and B blocks (a CB row) and a horizontal edge adjacent to a DA row. Similarly 
we will have a vertical edge adjacent to an AB column and one adjacent to a CD column. 

It is easy to convince ourselves that even with generic r b.c. we can repeat the same sequence 
of splitting and gluing, following the same "path" joining the 4 sub-lattices of Cl namely D — > 
C — > B — > A. In the bulk, namely where the sets D, C, B, A do not touch the boundary, we get 
the same results as in the case of periodic b.c. For the blocks close to dA p we get the following 
modifications: 

i) The various sets C,B,A of the bulk are substituted by their "truncations in A p " namely by 
C n A p , B n A p , An A p with the proper r b.c. on their part touching dA p and 0, like in the 
bulk, otherwise. 
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ii ) The various probability measures //^ , fi% k are defined similarly to what is done in the bulk 
with the difference that, in their definitions, the terms corresponding to partition functions 
on regions lying totally (resp. partially) outside A p are absent (resp. truncated); moreover 
the configuration on which they depend parametrically: a, (3 in fi^ k ; a in /i|j fc may contain 
r; notice that fJ-A kl stays unchanged. 
Hi) Some of the bonds, close to dA p , are consequently modified and their weights can depend on 
r. By an abuse of notation, we still denote them by , cj^\ \ A^\ 
Indeed the splitting operation is very similar in the bulk and close to the boundary; the true 
difference is the following. When we have some term produced by a splitting that, following the 
"bulk rule", we would like to glue with some other term outside A p or coming from A p , simply we 
omit the gluing and in this way we construct some new domains just consisting in the parts of the 
corresponding bulk domains (C,B,A), lying inside A p . 

Let us describe an example. Suppose that the upper horizontal side of dA p is adjacent from 
the exterior to a CB row (which, indeed, is the case with our choice of the location of Ap). After 
integrating over 5 variables and splitting like in (3.16) we do not glue on the blocks Ck 3 sitting on 
the top row like in (3.19) but we make an analogous operation combining the term Zo k3 (coming 
from the splitting on the Dj~ 3 block in A p ) with the self-interaction in Ck 3 and its interaction with 
the exterior configuration r. In other words we use a formula analogous to (3.19) but without the 
term Zo k +e2 which, now, is absent. In this way the set corresponding to Ct 3 in the bulk, just 
consists, now, of Ck 3 U Dk 3 - Accordingly we define fj,^ by omitting the factor Zr> k3+e2 (jk 3 ) i n its 
definition. When we continue with the splitting on the horizontal direction and the gluing, say, on 
Bk 3+ei we end up with the construction of a set, playing the role of J3fc 3+ei , obtained by removing 
from Bk 3 + ei the two D blocks exterior to A p where the "top" b.c. are given by r whereas the other 
b.c are still given by the reference configuration like in the bulk. Of course also the error terms 
(of <E> or \Jj type) are, accordingly, modified. 

In this way we can repeat the transformation of our system into a polymer gas. We just have 
to introduce the obvious modifications in the terms appearing in the expression of the partition 
function of the reference system n (see (3.55)) as well as in the bonds cj^\ B^\ 

A^p close to the boundary and in the measures , n% k when Ck 3 ,B k2 happen to be adjacent to 
the boundary dA p ; then, accordingly, we modify the definition of the polymers and of their activity, 
Cr = (r(Qr) (see (3.57)) which, now, will in general depend on the location of the polymer and on 

the b.c. r. Anyway if dg [r, I p ^j > d the activity of R is the same as in the bulk and does not 
depend on r. 

Proof of Theorem 3.1. Let us take the logarithm of (3.54). By using (3.55) (3.61) and (3.66) we 
get the following expression for the renormalized Hamiltonian. 

H^\n) :=log [Z^l] = J>g [z^ ((0))] - J>g [Z Fk ^ (0)] - £log [z^ ((0)) 

fci k 2 k 3 

k 

+ log [Z Dki ((0))] + E ^(Ri, ..,Rn)U C Rl (n 5i ) 



fc4 fc ^ 1 fli,...,fl*:fl i CA p i=1 
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We have 



H^\n) = const + £ f'^KH £ *S$'""(m*) 



(3.68) 



where, with G -4 P , -Bfc 2 G S p , Cfc 3 G C p , -Dfc 4 G V p and de(X,Ip) > <i (see the above discussion 
for dt(X,Ip) < d), we set 



(£,T),sr 



(mi) := 



and 



log 


"1 z< M < 


= mi,Q e (i) C A fel ) 


if X : U ieX QK0 


= A fcl 


-log 


A 2 ^ (Ml 


= m^Qtii) C F fc2 ) 


if X : U ieX <2^) 


= ^ 2 


-log 


Cfc 3 ,z 


= mi,Q e {i) C (7 fe3 ) 


if X : U ieA: Q*(i) 


= Ck 3 


log 




= mi,Q e {i) C D fc J 


if X : U iG xQ£W 


= D ki 


, 






otherwise 








fc 






ill,..., 









(3.69) 



(3.70) 



By the above discussion on the dependence of the activity on the boundary condition, for 
each X CC £e such that de(X,Ip) > d, ^x'^ ^ s independent of r. Therefore the limit in (3.6) 
exists and is actually reached for a finite p. Finally the estimate (3.7) is a direct consequence of 
(3.70) and Proposition 3.2. □ 
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4. The mult i canonical measure. 

Given a positive integer I and a volume A CC C of the form (2.10) we want to study the multi- 
canonical state which is obtained from the multi-grancanonical one by fixing the total number of 
particles in each cube Qe(i), i €. I. Let thus N_ = {Ni,i G /} be the random variables defined in 
(2.12) and, given n = {n^ = 0, • • • , \Qt\, i G I}, the multi-canonical state v\ n is given by 

"a,„(-) ■■=»l,A-\K = n) 

which, in the RG context, represents the constrained model. Note that v\ is independent on z. 

4-1- Thermodynamic relationships. 

We need to compare the multi-canonical and multi-grancanonical state. We start here by discussing 
some thermodynamic relationships between them. With respect to the usual treatment we work in 
finite volume and take advantage of the strong mixing condition to obtain explicit bounds. 

Let the volume A be of the form (2.10) for some I CC £e and fi\ z be a multi-grancanonical 
state satisfying Condition MUSM(.A). Introduce the map A 1 3 p T {z) G [0, l] 1 defined by 



\Qi 



-Ml,,{Ni), iei 



(4.1) 



Proposition 4.1. For each I CC £e and each closed C C A there is a constant C > such that 
for any boundary condition r, any z_£- C 1 and all I multiple of £ 



(4.2) 



1 + 



< Cz 



Q r e (i)riQ r e (j) 



-d(Q e (i),Q e (j))/C 



\Q 

pfU)-pJ(z) 



< n^i_ e -d(x,Q e (i))/c 
~ \Qe\ 



(4.3) 
(4.4) 



The proof of the lower bound in (4.2) is based on the following Gaussian bound on the 
characteristic function (see [DS4,§2.3] and [Y,§9]) which will be extensively used in the sequel. For 
t G R |J| , we use the notation (t,N) := Eie/^^i- 

Lemma 4.2. For each I CC £e there is a constant C > such that for any I and t G [— tt, 7r]l J l 



1 1, 



^A, £ (exp{i(t,iV)}) < exp<| -~\Qt\^Zi\l, 



CT 



(4.5) 



iei 



Proof. Before starting we stress that the proof is based only on the finite range and boundedness 
of the interaction and does not use Condition MUSM(«4). 
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Let A' C A be a subset of a sub-lattice of C with spacing parameter larger than the range r of 
the interaction. This means that for any x,y G A' we have d(x,y) > r but nonetheless |A'| > |A|/C 
for some constant C = C(r) > 1. If we set Q' e (j) := A' n Qe(j), we then have 



,i(t,N) 



r ( 171 

m*(do/&,, n n *** n 



< sup 

C 



a* In n 



n n sup|/,i, £ ( e ^K)) 



i=i^eQ^(j) 



since /U^, z is a product measure. 

Let p x (C) '■= P\i z (Vx = !)• Since the interaction is bounded we get, for some constant 
C = C(\\U\\) > independent on x, z and (, Zj/C < p x (() < Czj for x G Qe(j)- A simple 
computation on Bernoulli variables shows now that for |t| < tt, Xj G Qe(j) 



1 1 2 
<exp<J---^ 



the bound (4.5) follows. 

Proof of Proposition 4-1- We first note that 

in-Pi te) = 



□ 



Let 



dz 



r,{i) 



Zj\Qi 



(z) := z (Ni;Nj] 



the lower bound in (4.2) follows by noticing that Lemma 4.2 implies the quadratic form estimate 

i,j€l i€l 

To prove the upper bound in (4.2) and (4.3) we instead use Condition MUSM(„4) to get 

\vi,j\< Yl 2 W,AVx,Vy) < ^A,AVx,V x )+Cz iZj J2 e~ d{x ' v)/C 



xeQe(i)nQ e (j) 



and that for x G Qi(i), by the same argument as in Lemma 4.2, p\ z (rj x ,r] x ) < Czi. The proof of 
(4.4) is analogous and we omit it. □ 

Let p z be the infinite volume Gibbs state associated to the (translation invariant) interaction 
(z, U) satisfying Condition MUSM(„4). We introduce the (one dimensional) map A3 z ^ p(z) G 
[0, 1] by p(z) = fi z (rj x ) and denote by p z(p) the inverse map which is analytical as a consequence 
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of the strong mixing assumption. Let finally B C [0, 1] be defined by B := p(A) where A is as given 
in Condition MUSM(^); we note B = [0, 1] if A = [0, oo). 

Recall that the map z i— > p T '^(z) has been defined in (4.1). We need an inverse map p i— ► 
z^^\p) defined for all possible boundary condition r. When ,6 is a proper subset of [0, 1] we take 
i large enough and define it on a subset of B. By using strong mixing and Proposition 4.1 it is easy 
to deduce that for each closed C C A 

lim pi M (z) = p(zi), uniformly for r G Q, z G C J (4.7) 

and that for each closed set V C and any £ large enough (depending on V) we have 

P'cfV'WGA') 

T 

Finally, by (4.6), the Jacobian of the map z i— ► p T '^(z_) is not degenerate uniformly in r and Let 
P C £> be a closed set and I large enough; we can therefore define the inverse map on the set T> 1 , 
i.e. the map V 1 3 p ^ z_ T (p) = z_ T '^(p) such that 

P r (z r (P))=P 

for any p G 2? J , 

When B = [0, 1] we can instead define the inverse map for any t. Indeed we have 
lim pl'^\z) = 0, lim p['^{z) = 1, uniformly for r G CI, {zj G [0, oo), j / i} 

which, together with (4.6), implies 

£ ^)([0,oo) / ) = [0,l] / 

We prove below some estimate on the Jacobian of the map p i— > 2 r '^- > (/c); in order to describe 
them we need some more notation. Let {uJh, h = 0, ■ ■ ■ , k} be a path on the rescaled lattice Cg such 
that de(u h -i,Vh) = 1, h = 1,- ■ ■ ,k. We introduce q{u h -i,u h ) ■= \Q r e (u)h-i) n Ol(^)l/|Q^I- 

Proposition 4.3. For each £ Z + , J CC £< and each closed V C £> there is a constant C > 
such that for any r G f2, p G I) 7 , x G <9 r A and a/I I large enough 

h £ h A (£) £ c (4 - 8) 



sup sup i J q(uh_i, u>h) + >, W (4-9) 
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Moreover 



zf{p)-zl{p) 



< gft [ e -d(x,Q,(0)/c + gup sup gup TT q ^ h _ M ) + 

1^1 1 j:xed r Q e {j)l<k'<k ^ 



(4.10) 



Proof. Let 



be the Jacobian of the map z t— ► p T (z). We split in in its diagonal and off diagonal part; J = ID + . 
where 



and note that from (4.2), (4.3) it follows D > 1/C, ||A|| < CtT 1 . 

In order to prove the bounds (4.8), (4.9) we need to invert the Jacobian X We use the above 
splitting and Neumann series to get 



- 1 (1 + AD" 1 ) 1 = B" 1 |^(-1)' ! (AD- 1 )' 1 + R k+1 J 

\h=0 / 



where ||Rfe + i|| < Cl~( k+1 \ Since D is bounded from below and Ajj is exponentially small for 
de(i,j) > 1, (4.9) follows easily from (4.3). 

To prove (4.10) we note that, by definition of the map p i-> z_ T (p) we have 



pf (^(£))=P[(1 T (P)) 



% 6 J 



(4.11) 



By using the invertibilty (uniform in r £ and £) of z i— > p T '^{z) and (4.7), it is not diffucult 
to see that (4.11) implies that, for £ large enough, z r (p) and z r (p) are in the same connected 
component of z7* (P 7 ) U z r (V 1 ) . 

On the other hand, by Lagrange theorem 

pfU 2 )-pfU 1 ) = E^-4(D-[^-^}] 

where z_e A 1 if z 1 , z 2 are the same connected component of A 7 . Whence, by usig (4.11), 



zf ( £ ) - ^(p) = ^ (r(D)r. 1 • P j (^(p)) - pf u t (p)) 



and (4.10) follows from (4.4) and (4. J 



□ 
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4-2. Comparison of ensembles in finite volumes. 

We here discuss the equivalence of multi-grancanonical and multi-canonical ensembles. We shall 
work in finite volume with the aim of obtaining explicit bounds as a consequence of the strong 
mixing assumption. 

Let I CC £g, and A as in (2.10). We want to compare the measures fi\ z and v\ n where 
the activity z is chosen, depending on n, A and r, as (recall that the fuction p z_ T (p) as been 
defined above) z_ = z_ T (n/\Qe.\), i.e. so that p A z (N_) = n. We have the following result. Recall that 
B = P (A). 

Theorem 4.4. Assume p\ z satisfies Condition MUSM(A). Then for each closed V C B, each 
I CC £e and each local function f, there is a constant C depending on the constants in Condition 
MUSM(A), V, \I\, diam ||/||, such that for any b.c. r, any n £ V 1 and all i multiple of l 

the following bound holds 



u A,nf ~ MA,z/ 



<C^- . (4.12) 



The proof of this theorem is based on the DS complete analyticity conditions [DS1], [DS2], 
[DS3]. Although originally formulated for arbitrary volumes their theory carries over to our strong 
mixing for regular domains as already remarked. 

More precisely we need the following condition [DS3, Condition lb] which is equivalent to 
SM(4>)- There is a constant e > such that for all complex interactions $ in an e-neighborhood 
of <E>, i.e. 

$ G e ($) := {||$ - $|| < e} 

and all finite volumes A as in (2.2) the analytic functions Z T K {<&) are non-vanishing. Moreover, 
there is another constant A 1 < oo such that for all <E>i, $2 £ C e (^ > ) we have the bound 



sup 



PX($i) - Pa($ 2 )| < A' \T n supp ($1 - $ 2 ) I (4.13) 
where the pressure P is defined by 

PX($):=logZX($) (4.14) 

and supp (<J>) := Ua : $ A /o A. 

Proof of Theorem 4-4- Since the b.c. r is kept fixed we drop it from the notation. We also assume, 
without loss of generality, that ||/|| is small enough. 

Step 1. We express here the difference between multi-grancanonical and multi-canonical states by 
introducing the Fourier transform of the indicator l/v= n . 

By definition of the multi-canonical state u^ n , we have 

m „ ( n _ ((/ ~ /*A,z (/)) %=n) _ MA, £ ((1 + u)1n=u) , 
"A,n (/) - MA,z (/J = TTT x - TT-j x 1 (4.15) 

MA,zQV; = n) MA,z(iV = n) 
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where we introduced u := f — ha,x (/) which has the same support as / and is mean zero w.r.t. 

We next introduce the perturbed probability measure dji\ z := (1 + u)dfiA, z - We regard it as 
the Gibbs measure w.r.t. an interaction <J> M . Since / is a local function, we have that <J> U has range 
bounded by max{r, diam (supp(/))}. Moreover, by taking ||/|| small (depending on e) we have that 
$ M G e ($). 

By taking the Fourier transform on the r.h.s. of (4.15), we have (recall that ^a j£ (N_) = n by 
the choice of z) 

(4.16) 

f d* e ^A(t,z)-i(t,AtA, £ iV) r e VA(*.£)-V>A(t,z) _ I] 

J|t|<7T - L J 

J|I|<7T - 

where, indicating with a superscript the dependence on the perturbation u and inside the paren- 
theses the dependence on the complex activity, we introduced 

Mt,z) := log/i A , £ (e^>) = P A ({ Zj e^} j€l ) - P A (z) (4.17) 

where the second identity holds by expressing the l.h.s. in terms of ratio of partition functions. The 
definition of i)\(t,z) is analogous, it is enough to consider the pressure of the perturbed interaction. 

Step 2. Here we estimate from below the denominator on the r.h.s. of (4.16). 
Let us introduce the variances 

v1 = vf\zf := nl^N.-N,) 

and note that from Proposition 4.1 we have C~ l Zi\Qg\ < vf < C Zi \Qi\. This bound will be used 
extensively in the sequel. 

We shall prove the following bound. There is a constant C independent on r, t and z such 
that for I large enough 

»A, Z (K = n) = 7 -l Jr / dte*^-^"^ > i TT ^— (4.18) 

where we recall z has been chosen so that //a j£ (TV) = n. 
By a change of variables we get 

Ma z (N. = Ik) = tt — ~o f dse^^'^-^'^^) 

where we used the notation s/v to denote the variables {sj/vj, j G /} 
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Let K be a large constant. We take advantage of the Gaussian bound in Lemma 4.2 to get 



L 



d s gV 1 A (s/v ,z) - i (sjv , fJ, A ,zN_) 



3j:KA(nVj)<\sj \<7TVj 



< 



L 



ds exp 



3j:\sj\>KA(-rTVj) 

< Ce- R2 / C 



1 1 

2C 



iei 



(4.19) 

By the above bound we can restrict ourselves to bounded s. We need however to treat 
separately the Gaussian scaling in which Vi diverges with I and the very low density case in which 
it remains bounded. Let M be another large constant (1 <C K <C M <C I) and introduce I g := 
{i £ I : vf > M}, I p := I \ I g . Let also s g := {si,i £ I g } (resp. s p := {si,i £ I p }); we 
use an analogous notation for z. We shall prove the following expansion on the logarithm of the 
characteristic function. 

tp A (s/v,z) - i {s/v,fi A ,zK) 

- £ Lisi/v, _ x _ iSj/ \ _ 1 J- (N 3 ;N f ) ^+Ra (s,z) (420) 



where 



.„ , / K 3 M 2 KM K 2 M\ 

sup \R A (s,z)\<C —= + — — + , — . + — — 

\s\<k \VM \Q t \ t/WI \Qt\ J 

Note that on the r.h.s. of (4.20) the first term corresponds to a Poisson limit for Nj, j € I p and to 
a (joint) Gaussian limit for Nj, j G I g . 

Postponing the proof of (4.20), let us first show that, together with (4.19), it implies the 
bound (4.18). It is enough to notice that if Z is a Poisson r.v. with mean A £ Z + we have 



-j- [ ds e {^-^s/u)X = u p rob {z = A) 

27F J\s\<WU 



u- 



A! 



By using the bounds vf > Zi\Q(\/C, fi A ,z^i < Czi\Q(\, Stirling's formula and estimating the 
Gaussian integral (recall (4.6)) we thus get 



ds exp < 



\si\<KA(-KVi) 



> 



C 



and (4.18) follows since we can make the remainder as small as we want. 

In order to prove (4.20) let us first expand in power series of s g and get 







V>A {s/v,z) =1p A (0,S p /v,Z_) + E ^T^A (0,S p /v,l) — 



iei„ 



Si 
Vi 



+ 2 2^ duduM Q ^^ v^ + R ^ l] 

i,i'€l g 



(4.21) 
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We note that by Condition [DS3,Ic], still equivalent to SM(£ ), 



^U,l) = \v T A,z,t_ (Ni)\ < C Zi \Q e \ (4.22) 
here \i\ z t denotes the complex measure defined by 

We remark that in [DS3,Ic] does not include Zi on the r.h.s. of (4.22). However, by the remark 
following Condition USM(^ ), we can easily verify that (4.22) holds. 

Recall that the pressure Pa(z) is holomorphic in an e-neighborhood of z. Therefore (see 
(4.17)) ipA(t,z) is holomorphic in a neighborhood of t = 0. By taking K/y/M small enough we can 
thus use Cauchy integral formula and bound the third order derivatives (w.r.t. to t) of tp(t,z) in 
terms of the first one. By applying (4.22) we get 

Si Sj Sfz | 



sup \R\(s,z)\<C sup mm{zi, z~, z k }\QA—^ — < CK 3 



\s\<K 



ViVjVk 



We next expand the other terms on the r.h.s. of (4.21) in power series of z p . Note in fact 
that for i G I p we have zi < CM/\Qg\. Let us consider the first one. We get 

d 

V'A (0,s p /v,z) = IpA (0,s p /v,z g ,Q) + ^2 Zj-Q^-^A {0,s p /v,z g ,Q) +R\ (s,z) 

Noticing that \ipA(t,z)\ < C\Qe\ and using again the Cauchy integral formula, we can bound the 
remainder as follows 

M 2 



\Rl(s,z)\ < C z J z f\Qe\ < c 



\Qt 



We next observe that V'A (0, s p /v,z g , 0) = 0. On the other hand, by (4.17) 

(0,s p /v,z g ,0) = (e^ - l) A Pa ( l3 , ) 



By the analyticity of the pressure (see (4.22)) we also have 

< C\Q t \ £ * < CM 



Pa{z 9 ,0) - —P A (z g ,z p ) 



dz 



iei v 



Since Zj-§-P A (z) = ha,z w e thus get 



^A (0,s p /v,z) = J2[e 



l)»A,ANj) + Rl (s, z), \R 3 A (s, z)\<C 
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We expand similary the other two terms in (4.21). For the second one we have 



7^A (0,s p /v,z) = — 4> A (0,s p /v,z g ,0) +R\ tt (s,z) 



dt 



where, by using again (4.22) and the analyticity of V>A (t,z), we have 

\Rii U,z)\ <C Zi \Q e \Y^ Zj < CMzi 

Furthermore, since by setting Zi = ipA becomes independent of Sj, 



d d 
t^a i°^ p /v,z g ,0) - — ip A (0,0, z g ,z p ) 



dti 



d 



q^^a (0,s p /v,z g ,0) -inA,z(Nj) 



< CMzj 



so that 



d 



jei g °' < v i r /, Vj 



sup \Rl(s,z)\<C^L= 



By the same argument we finally have 



d 2 d 2 

^A{0,s p /v,z g ,z p ) = — — ^ K {0,s p /v,z g ,0) +R Ati j(s,z), |# A ,i,i'U>*)| < CMziAzi 



dti&ti' 
Moreover, as before, 



(4.23) 



^—^A (0,s p /v,z g ,0) - ^—^a (0,0,z g ,z p ) 



d 2 



dtjdtji 



iPa {0,s p /v,z g ,0)+ t x A ,z(N i ;N i ,) 



< CMzi A z v 



which gives us 



i'ei g 



Vi Vi 



where 

K 2 M 

sup \R\{s,z)\<C-= 

\§\<k v\Qt\ 

The proof of (4.20) is now complete. 

Step 3. We finally here estimate from above the numerator on the r.h.s. of (4.16). 
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Let Ke := log \Qe\- We make the change of variables t = s/v and use Lemma 4.2 (which holds 
also for the perturbed measure n\ z ) to get 



d s e '4>A{s/v,z)-i(s/v, t j, At£ N_) e 4>%(s/v,z)-ip A (s_/v,z) _2 



3j:Kif\(iTVj )<| Sj \ <^Vj 



<-j 

J3j:K, 



ds 



3j:KiA(-7TVj)<\sj | <nvj 



+ 



< Ce- K "/ C < C 



\Qt 



We can thus consider the case | Sj \ < Ke A (irvj ) . 

Since Zi\Qi\/C < vf < Czi\Qe\, either Si/vi or Zi is small. We can therefore apply the bound 
(4.13). We get 



P«({ e ^/«^.})-P A ({e^/^^.}) 



< C 



We next expand the difference ^~a{s/v,z) — iPa(s/v,z) in power series of s. Since [i\ z (N^) 
VA,z (N k ) = MA,z (/; N k ), we get 

^aU/v,z) - iPa(s/v,z) = i^2^A,z(f;N k ) — +R{(s,z) 

kei Vk 



where 



i,jei 3 



&l(t,z)-M±,z)} 



t_=l/v V i V 3 



We note that, by (4.17), 



J- [m_,z) - Mt,z)} = iz k e*> ^ [PXW) - Pa(z')} 



By the analyticity of P]^(z') — Pa(jl'), for t = s/v we can bound the r.h.s. above by Czk- We thus 
have 



\Ra(s,z)\ <CV ^Az,^ < C 

1 ' 11:11- 



ViVj ' \Q e \ 



As \nA,z{f'-,Nj)\ < Czj, for |s| < Kg we finally have 

ew{^l(s/v,z)-^ A (s/v,z)}-l = i ^/i A , £ (/;iV fc )-+ J RiU,l), |«i(s,l)| <C 



fee/ 



(4.24) 



By Lemma 4.2, we have 



ds e ^(»/«,5)-i<»/«,MA > ,(JV))_ R 2^^ 



IsjI^^aCttUj) 



< c 
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To conclude the proof we consider separately each of the other terms on the r.h.s of (4.24). 
We want to show that, with a small error, the function tp(s/v,z) — i (s/v,[ia,z (N.)} is even in Sk', 
hence the integral vanishes by symmetry. We thus expand ip(s/v,z) as follows 



1 d 2 s s 

tp(s/v,z) - i {s/v,fi A ,z (AT)) = - ^ ^ tp(s/v,z)^"- : 



2 ^— ' dtjdtji 

3,3 



by letting s^ := {s i} i €. I\ {k}}, we have 



d 2 



dtjdtji 



Hl/v,z) = -B^d^) + R* M, z), B^,(s^) : 



dtjdtji 



tp(0,l {k) /v,z) 



and, by (4.22) and the analyticity of tpA, 



\Rl,j,fU,z)\ < Czj Azj'lQtl^- 



Whence 



^{s/v,z)-i(s/v,n^{N)) = -\ £ \RU(s,z)\<C 



3,3'ei 



We next use the bound 



Vk 

(4.25) 



e «l fc (£.£) _ i| < (i + L<*(»,*)h \R\ k (s,z)\ 



and (4.25) to get 



\sj\<K e A{TVVj) ' v k 



J 

J\s. 



< I ds 

'\sj\<K t A(Trvj) 



+ 



|#A,fcU,l)| 



(4.26) 



where, recalling that \ha,z (/; A^) | < Czfc and u| > Zfc|Q^|/C, 



Rlk(s,z) := MAi£ (/; AT fe ) ^< fe U,l), « fe U,l)| < c|| 



By applying again Lemma 4.2 we have 



1, 



ds 



\sj\<K e A(TVVj) 



3 Va (s/v,z) 



\R 5 A,k(s,z)\<C 



\Ql 
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It now remains only to estimate the other term on the r.h.s. of (4.26). Let Mg := £ 1 / 4 and 
introduce L [ g k) :={»€/ \ {k} : vf > M £ }, I { p k) := I \ ({k} U . We have 



B 3 ,Ai {k) ) = - 



d 2 



dtjdtji 



msP/v,z) + Rl k ^,(s,z) 



where 



SUP \R 6 A:k>hJ ,(s,z)\ < CZj A Zj'\Q(\ ~7~ — Cz j A z i'-rfr 

\s\<K e :^ri,_\ Ul y J-yJ-t 



iei<-k) a 



K f 



so that, by using also (4.23), 



E s »'d (fc) )^7= E ^MM^ + Rlk&z) 



where 



sup \R 7 A , k (s,z)\<c( 

s\<K e \ 



Kf + KjM, 



fMt \Q t \ 



Hence, recalling (4.6), 



L 



ds 



\s j \<K e /\(lTV j ) 



which concludes the proof. 



_IV B ,(c (k) ) SiSj ' 



K fe u,i)|<c 



□ 



^.3. ZocaZ Central Limit Theorem with multiplicative error 

In order to obtain the convergence of the short range part of the renormalized potential to the one 
of independent harmonic oscillators we need a local central limit theorem which will allow us to 
compute the asymptotic behaviour (as £ — ► oo) of the r.h.s. of (3.69). Since we are interested in 
the logarithm of the partition function we do need a local CLT in which the error appears in a 
multiplicative way. It can be proven by applying the theory of moderate deviations as developed in 
[DS4]; alhough these results are stated only for very high temperature, the proof is based only on 
the analyticity properties of the thermodynamic functions which hold under Condition MUSM(„4). 

Let us recall that fi A z is the multi-grancanonical state in a volume A C C C of the form 
(2.10). We denote by v ^ = v T <( £ \z) the covariance matrix of the total number of particles in each 
cube Qe(i), i.e. v T ^ e \z)ij := \i\ (TV- Nj), where Ni has been defined in (2.12). We have the 
following local central limit theorem. 

Theorem 4.5. Let U satisfy MUSM(A) and pW = p T '^(z) := fi T A ^(N) /\Q e \. For each A of the 
form (2.10) and z G A, e > there are constants 5 = S(z_, L, e) > 0, C = C(z_, L,e) < oo such that 
for any integer £ we have 



VA,AK = n) = [(2tt)I 7 I detvW 



"expj-l^n-pWlQ,!),^))" 1 ^-^)^!) 
x{l + R A (n)} 



(4.27) 
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where 



sup 



sup 

n : 

V f >IQ^l|<|cW 2 / 3 - E 



\Rl(n)\ < C m 



(4.28) 



This Theorem is essentially contained in [DS4]; however to make the paper selfcontained we 
give below a brief sketch of the proof. Given n we let C = C M \n) be defined by C := z T ' (e) (n/\Qt\) 
where we recall the function p ^ z T '^(p) has been defined in Section 4.1. We also recall the 
pressure has been defined in (4.14). We have the following Lemma. 

Lemma 4.6. Under the same hypotheses of the previous Theorem, there are constants = 
£ od;-0 > 0, C = C(z,I,Eo) < co such that 



where 



and 



MA,zQV = n)= [(27r)l J ldet^(C)] * exp {-/ A (n)} (l + Rl(n)) 
Ia(u) = Il,An) ■= 5>ilog^ - [PKO ~ Pl{z)] 



(4.29) 



iei 



sup 



sup 



Rl(n) 



< C 



1 



\Qi 



(4.30) 



™-p (£) IQ^I <enl<3£l 

Sketch of the proof. By definition of the multi-grancanonical state p\ z we have 



MX,,Qv = n) = n \r 



z T (Q 



-•fA(n). 



1 



(2vr)l 7 l 



l*l<» 



If we take e small enough, |n — p^|Q^|| < implies that (£, i7) satisfy SM(£ ) for some 4> = 

(l) £ o ) • In order to conclude the proof it is then enough to make the change of variables t j = 

s i/ \f^fl-> use Lemma 4.2 to estimate the tail and expand logp\ ^ (e^-'— >) up to the third order, 
using analyticity to estimate the remainder (see Section 4.2 for analogous computations). Note in 
fact that, by the definition of ( we have p\ ^ (TV) = n. □ 

Sketch of the proof of Theorem 4-5. By applying Proposition 4.1 we have 



sup 

Ten 



sup 

— '■ 

U-pW|Q^l|<IQ^I 2 / 3 - E 



w (/) U)-« w (C) <C|Q/ 



|2/3-e 



which, together with the bound (4.6), implies 

(det^C))" 1 = (det V «>(*))~* (l + i? A ' (1) (n)) 

i? A ' (1) (n)<C7 1 



where 



sup 

Ten 



sup 

n : 



IQd 1/3 
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On the other hand, by the analyticity (uniform in £) of the thermoduynamic functions, we 
have (see [DS4 Eq. 1.2.15]) 



Ia(u) = U(n-p w \Qt\) , (vU®)' 1 (n-p w \Qt\)) (l + R? 2) (nj) 



where 



sup 



sup 



< C 



\n-p ( -^\Qe\\<\Qe\ 2 / 3 - e 

in which we have used again that ^ (TV) = n. 



\Q< 



\3e 



□ 



5. Gibbsianess and convergence 

In this section we conclude the proof of the main results. First, by applying the comparison of 
ensembles, we show the constrained models satisfy a finite size effective condition uniformly in the 
constraints. Secondly, by applying the local central limit theorem, we prove the short range part of 
the renormalized potential converges to the potential of independent harmonic oscillators. Finally, 
when the global condition GMUSM holds, we verify that the renormalized measure (defined 
directly in infinite volume) is Gibbs w.r.t. the potential constructed in Section 3 (obtained via a 
thermodynamic limit). 

5.1. Finite size condition for the constrained models 

We consider the BAT obtained by partitioning the original lattice C into cubes of side £, C = 
Uiec e Qe(i)- Let fj, z be the (infinite volume) Gibbs state of the original system at activity z. We 
then introduce the constrained system by fixing the total number of particles in each cube; it is 
described by the conditional (multi-canonical) measure we introduced in the previous section. 

We want to show that, provided Condition MUSM(„4) is satisfied, the local specification 
associated to the multi-canonical state v T K n satisfies (3.5) with 5(£) = C/£. We shall consider £ 

to be an integer multiple of £ . Recall that B = p(A), L = d£ and vf> = (\Qe\V) A n (see 
Theorem 3.1). 

Proposition 5.1. Assume the interaction U satisfies MUSM(A). Then for each closed setV C B 
there is a constant C such that for all L the following bound holds. 



sup sup sup sup sup 

ieC e k=l,-,d AeP r (fc) (i) n€V (e) CT ><> T 

Li K ' A 



Za~ 



s ((?(*■+>, *(*■-), t)Z A| „ (C (fc ' +) ,C (fc '- ) ,r) 
i2 (<T(*.+), C (fc,-), r )Z A ,„ (C( fc .+>,<7( fc .-\r) 



< C j (5.D 



Lemma 5.2. In the same setting and notation of the above theorem, there is a constant C such 
that for any A C A for which d (A, d^ k, ~^A) < r, diam(A) < r 

sup sup sup sup sup Var [u A A ' ,v\ A ' < C— — (5.2) 
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Postponing the proof of the Lemma, we show how it implies the main estimate. 
Proof of Proposition 5.1. Let us first show that (5.2) implies the following condition 



sup sup sup sup sup sup 

ied k=i,-,d xeP (k) d) xed^-)A ne v {e) CT 'C> T 



We have in fact 



Z A ,n {a^+\r^r) Z A , R (C (fc ' +) ,r,r) 



Z A , R (C( fc ' + ),r-,r)Z A> „ (a( fc '+),r,r) 



< C 



1 



\Qt\ 

(5.3) 



Z A ,n (C (fc ' +) ,T*,r)Z A ,„ (a( fc '+),r,r) 



- 1 



Z K n (a (fc ' +) ,r",r) 
Za, r (a( fc .+),r,r) 

Z A; n (a (fc ' +) ,r",r) 
Z A ,n (a( fc .+),r,r) 



^A,n(C (fe ' +) ,r,r) Z A) „(a( fc '+),r,r) 



^A,n (C (fc ' + ) , T* , T) Z A> „ (<T (*. + ) , T* , T) 



' T (^)-< >T (^) 



(5.4) 



where 



K(rj) ■■ 



— p -[Ha(v°AT x )-H a (v°At)] 



is a local function with support contained in an r neighborhood of x. Since the first factor on the 
r.h.s. of (5.4) is bounded uniformly and the same holds for \\hl\\, (5.3) follows from (5.2). 

An easy telescopic argument shows (5.3) implies (5.1). Indeed, for any two configurations 
(( k '~\ o~( k '~\ differing only on d^ k '~^A, we can find a path {?7;}z=o,---,M of length M < r ■ (3L) d ~ 1 
such that 770 = o~( k '~\ rjM = £( fc '~) and r\i differs from r]i-\ at most in one single site x € d^ k '~^A. 
We then write 

Z A ,n (a^ + \a^-\r) Z A ,„ (_^\ r) 

Z A ,n (<7(*.+>, C(*>->, r) Z a ,„ (C( fc '+), r) 

_ TT Z A ,n (<t<*- + U-i,t) Z K& (C {k ' + \m,r) 
l\ Z A ,„ (a( k '+), m ,r) Z A , R (C( fc '+),^_i,r) 



and use (5.3) to get (5.1). 



□ 



Proof of Lemma 5.2. Let us recall that Var (fi, v) = sup||j|| =1 \fif — vf\. Let / be a local function 
with support contained in A. By Theorem 4.4 we have 



< c 



+ 



V A ,z}f V A ,z 2 f 



where z_ a = z a (A,n, ( a ), a = 1, 2 is chosen so that ^° z0 QV) = n. Since Ci differs from ("2 only on 
9(fc,+)A, by Condition MUSM(^) we now have 



C2 



< Ce-< A ' d(k ' +>A )/ c 
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On the other hand, by Lagrange theorem, for a suitable z, 



ieA 



By the exponential decay of correlations we have 



a*L (/; * 



< Cz,e- d ^ A ' Q ^/ c 



the bound (5.2) is thus obtained by applying Proposition 4.3 to estimate \zf — z}\. Note in fact 
that d(A,d^A) > d£-r. □ 

5.2. Short range renormalized potential. 

In this section we consider the limit i — > oo of the short range part of the renormalized potential. 
By applying Theorem 4.5, we prove the necessary estimates. This would also allow us to conclude 
the proof of Theorem 2.3. 

Proposition 5.3. Recall that the short range part of the renormalized potential $> x have been 
defined in (3.69). We introduce 



(5.5) 



where 

'+lif X = A kl ,D ki 
S(X) := | -lif X = C k3 ,F k2 
otherwise 

v 

Then the renormalized Hamiltonian can be written as 



Moreover there is a constant a > such that 



lim sup $j' sr (mx) 



for any X CC £/, |X| > 2 



(5.6) 



(5.7) 



(5.8) 



Note that Theorem 2.3 follows directly from Theorem 3.1 and Propositions 5.1 and 5.3. 
Proof of Proposition 5.3. By using (4.27), for each V = Ui e xQe(i), recalling that rrii = (nj 



p\Qe\)/Vx\Qe\, w e have 

log Hy z (Mi =mi, ie X) =const - - ^ 



+ log [l + 4(mx)] 



(5.9) 
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Therefore, by (3.69) (where the boundary condition is r = and de(X, I£) > d), we have 



= S(X) U(n-p«>\Q t \) , (JOY 1 (n- p«>\Q t \)) - \ £ 



2 \ V- ^ '"'7 'V ) \~ '- '"'7 / 2 



(5.10) 



+ io g + 



Indeed, to get (5.7), it is sufficient to observe that 

i) given a block Qe(i) contained in A p , the corresponding one-body renormalized interaction 

— \m1 appears only in one term $f' ,sr (mj() with X = for one and only one G A p 
with S(X) = +1 

ii) given a block Qe(i) contained in B p , the corresponding one-body renormalized interaction 

— \rn1 appears in two terms $j' sr (mj() with X = Ak 1 with S(X) = +1 and in one term 
with X = F k2 with B k2 G B p and = -1 

Hi) given a block Qe(i) contained in C p , the corresponding one-body renormalized interaction 

— \m\ appears in four terms $f' sr (mi) with X = A kl with S(X) = +1, in two terms with 
X = Fk 2 with Bk 2 G £> p and = —1 and in one term X = Ck 3 with Ck 3 G C p and 
S(X) = -1 

iv) given a block Qi(i) contained in V p , the corresponding one-body renormalized interaction 

— \m1 appears in two terms $f' sr (mj() with X = with S(X) = +1, in two terms with 
X = Ck 3 with Cfc 3 G C p and 5(X) = — 1 and in one term X = Dk A with Dk A G £> P and 
5(X) = +1 

Performing the different cancellations in the four sub-lattices A p , B P ,C P ,V P we easily get (5.7). 
Finally, to prove (5.8), we note that by Proposition 4.3 we have 



\Q e \U e A 1 -5 i , 3 

v / t.j 



and, by strong mixing, 



< cr 1 



Hence the bound (5.8) follows from (5.10) and Theorem 4.5. □ 

5.3. Gibbsianess of renormalized potential 

We show here that, provided Condition GMUSM holds and i is large enough, the renormalized 
measure is Gibbsian w.r.t. the potential which has been constructed in Section 3. We 
have in fact the following result. 

Proposition 5.4. Assume Condition GMUSM holds and define the renormalized potential 
as in Section 3. Then the renormalized measure /JP is Gibbsian w.r.t. i.e. 

(m/| m 7 c) = ^ i T , ii%> a.s. (5.11) 
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Note that Theorem 2.2 follows directly from Theorem 3.1 and Propositions 5.1, 5.3 and 5.4. 
Indeed GMUSM implies B = p([0, oo)) = [0, 1] 

Proof of Proposition 5.4- We recall the random variables Mj = Mj(r/) have been defined in (1.3). We 
introduce the two families of cr-algebras: J-\ := cr{ri x , x £ A}, Ac£, and := a{Mi, i £ /}, 
I C C t . For I CC Ct and F : i-> R let us first prove that 

/4*> (F(m 7 )| ^) = ^ (f(Mj)\ Ff)) , M W a.s. (5.12) 
let G be a local function measurable w.r.t. ^f); by definition of the measure we have 



fi 2 (F(Mj)G(M Ia )) = ^) (F(mj)G(m IC )) = J d^(m) G(mj<>W (F(m 7 )| ^ 
on the other hand, 

M*(F(Mj)G(Mjc)) = fduM G(M lB (ri))fi z (f (Mj( v ))\ 7$) 

r d»<P(m) G(m 7 c) M , (^(MKr?))!^) 



which proves (5.12). 

Let V = ^l i€ yQe(i) CC £; we note that for I C V" we have 



= Hz (fi z (Mj = mi \ 4^ V | ^£>) 



(5.13) 



on the other hand, by definition of the renormalized Hamiltonian and the corresponding potential, 
see Section 3 



>. xcv x " v v\i I ( 
»v,z[Mi = m I \M vv = m vV ) = ^ ^ ) —- Y (5.14) 



exp < 






D roie nw ex P < 





Since Condition GMUSM holds, by Proposition 5.1, (3.5) is satisfied with V = [0, 1] and therefore, 
by Theorem 3.1, the r.h.s. of (5.14) converges, as V | £, to the r.h.s. of (5.11) uniformly in r and 
m. By using also (5.13) and (5.12) we thus conclude the proof. □ 

A.l Proof of USM(„4) => MUSM(i) in dimension 2. 

Let RL,3L(i) be the rectangle with vertical and horizontal sides L, 3L, respectively, and which 
is centered at 

The fact that we only consider this rectangle with longer horizontal side does not represent, 
of course, a loss of generality and is made only to fix notation. 
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For M an even integer, M/2 and L odd integers, we write: 



L = ML ; R iXL = R LfiL ((**-l, 

again the choice of the center is made to fix notation and does not constitute a loss of generality. 
Recall that since M is even and L is odd the center of Qiii^f 1 , ^P 1 )) is in 4f). 

We set Ri t 3i = Q\ U Q\ U Q r L where by Q^, Q|, Q r L we denote the left, central and right 
L x L squares, respectively, contained in Rl,sl- 

Consider a 2D lattice gas with an interaction satisfying USM (A) for some A C [0, oo). We 
start noticing that from the validity of USM(„4) it is immediate to deduce that for each z G A 
there exists an integer Lq such that the following condition 



sup sup 

<j,t€Q ie{l,2} 



Z v (cr( i '+),r( i -),r) Z v (r( i '+),a( i -),r) 



< e(2) (.41.1) 



is verified for V = Ql {})i RL ,3L (i) m the homogeneous activity case. This, together with the 
results of [O], [OP] establishes the equivalence of USM and CI in the homogeneous activity case; 
this result is valid in any dimension. 

Now, given a closed set C C A suppose that we are able to prove the existence of L such 
that: for all z, z' G C, if we consider our lattice gas enclosed in V = Ri^L with activity z' in Q l L 
and z in Q C L U Q r L (i.e. we take the same activity both in Q| and Q r j), then, calling Zv, z , z '{t) the 
corresponding partition function with r boundary condition, we have: 



sup sup sup 

cr,T ie{l,2} y£d(i, + )v 



Zv,z,z'(cr y ,ay' , r )Z v ,z,z'(T y , T y >, r 



Zv,z,z'(<Ty, T y' ^) Z V,z,z'(Ty, <V , 1") 



with e(Z) going to as L goes to infinity; then, using methods and results of [O], [OP] it is easy 
to get MUSM(„4). Indeed in the two-dimensional, multi-grancanonical case, to get strong mixing 
condition using effectiveness of some finite-size conditions for volumes of the form (2.10) with £ 
sufficiently large, it is sufficient to verify: 

(i) (Al.2) for V = Qi(i) and V = Ri 3 i(i) with uniform activity in V arbitrarily chosen in C, 
and 

(ii) (Al.2) for V = Ri y si{i) and activity z' in Q 1 ^ and z in Q c ^ U Q£ uniformly for z, z' in C. 

In the homogeneous case (i), as we noticed before, if, given C, L is the size for which SM(L ) 
holds uniformly in C, as prescribed by USM (.A), then, for L sufficiently large (Al.2) holds for 
V = Qz{i) and V = Ri^i{i) for each (constant in V) activity z G C. Then 2.1 will follow from 
next Proposition Al.l 

Proposition Al.l. Suppose that Condition C1^ 2 \V) holds for any V = QL (i), RL ,3L (i) con- 
tained in one of the three squares Q l L , Q C L or Q r L ; then, for M = sufficiently large, (Al.2) holds 
forV = R LM . ^ ' 

Proof. 
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We make a geometrical construction similar to the one introduced in [0], [OP] and used in 
Section 3 to compute, via cluster expansion, the renormalized potential. We recall that we denote 
by C our original lattice Z 2 whereas we denote by C l the L -rescaled lattice: we partition L into 
cubes of side L . We write: 

From now on we will mainly consider the Lo-rescaled lattice; our unit length will be Lq. In other 
words we will use the distance cIl . The "bricks" of our construction will be the blocks Ql 
or i?L 0i 3L n and the original length-scale will enter only when considering some properties of the 
partition functions in the regions Ql or Rl q ,3L that we use as input of our perturbative theory. 

Let ei , e 2 denote, respectively, the horizontal and vertical lattice unit vectors in Cl ■ &i = 
(1, 0), e2 = (0, 1). Following definitions and notation of Section 3 we further partition Cl into four 
sub-lattices: 



where: 



£-2L„ '■= {i = (h,i2) e £ Lo ■ i\ = 2jx,x 2 = 2j 2 ,for some integers 2/1,2/2} 



B 

^2Ln 



£L a + ei + e 2 = C$ Lo + e 2 

£2L„ + ei = + e 2 = £f Ln + ei + e 2 



(A1.3) 



We also set, for i G £l : 

:= Q Ln (2i), Si :=Q Lo (2i + e 2 ) C, := Ql„ (2i + ei + e 2 ) D t := Q Lo (2i + ei ). (A1.4) 
Then we can partition V = Ri 3 i into the union of the L -blocks of the four types: A, S, C, D: 

V = A V UB V UC V UV V 

where 

.Ay :={A i :i = (i 1 ,i 2 ) G £ Lo : N < (M/2-l)/2, |n| < (3M/2 - l)/2} 

and similarly for By , Cy , 2?y . 

We have that the left block on the bottom is an A-block whereas the right one on the top is 
a C-block. 

We denote by a, a generic spin configuration in Ai. ai G { — 1,+1} L °. Similarly for (3i,ji,6i. 
We simply denote by a, (3, 7, S the configurations in Ay , By , Cy , T>y , respectively. 

Notice that we have used the same notation (with a very similar meaning) as the one we used 
in Section 3 to describe "multi-canonical" block variables. 

Consider the "column" V\ namely the rectangle with basis L and height L placed at the 
left-hand of Q| , adjacent, from the exterior, to Q 1 ^: 

. L Ln + 1 L Ln + 1 LLn + 1 L L - 1 , 

V l = {(x 1 ,x 2 )€C:-- + -^- <x 1 <-- + ^ T - + L , -- + ^— <x 2 <+- + ^— } 
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we decompose Vj as disjoint union of A and B blocks: 

V l = Ai U Bi 

where 

Ai:=A v nv t , B r .= B v nVi 

We have: 

Ai = U ieI i A Ai 

where 

I l A = {{hM) : n = -(M/2-l)/2, |i 2 |<(_-l)/2} 

similarly for £?; 

We write .Ay = Al) Ai, By = B U Bf, in other words A, B denote the union of A and B 
blocks, respectively, which belong to V = Rz^l but not to Vj. 

We will repeat almost the same computation that we made , in the multi-canonical framework, 
to compute the renormalized potential. Namely we adopt the same strategy based on a block 
decimation procedure over the sequence of sub-lattices D, C, B, A. 

The main difference here is that we will treat in a different manner the region in V = Rt,zL 
adjacent to the boundary between Q C L and Q l L . Here we will exploit the fact that this boundary is 
one-dimensional. 

Indeed we will see that the system of the surviving a-variables in Ai, after decimation on 
5, 7, (3 and a in A, gets an effective interaction which is exponentially decaying with the distance 
and uniformly bounded in norm. The resulting one-dimensional system, regarded on a sufficiently 
large scale, is in the weak coupling region and from this it easily follows a weak coupling between 
opposite horizontal sides of V = Ri 3 i so that condition C3 with an infinitesimal £3 is satisfied for 
V. 

We want to perturbatively treat, similarly to what we did in Section 3, the partition function: 

Z^:= WWII))) 

where 

A:AnVV0 

and we recall that we are using the notation: 

V := Rl,zli t £ = boundary condition outside V 

Given r, r' G £lv c an d x, y belonging to the set of conditioning sites above the upper side 
and below the lower side of V, respectively, we want to consider the ratio 

Z V:Z:Z ,(t^t^t)Z V:Z:Z ,(t) 
Z V:Z:ZI (t^t)Z V:Z:ZI (t^,t) 
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(.41.5) 



where (r,T x ,Ty), (r, t' x ) ,(t, t' v ) are the configurations obtained from r by substituting r with t' in 
{x} and {y}, respectively. 
The perturbative expression that we will obtain for Zy will show an almost factorized depen- 
dence on boundary conditions in opposite horizontal faces so that we will be able to show that the 
quantity 



Z V}Z}Z >(T x ,T)Z VyZyZ ,(T',T) 



(A1.6) 



can be made arbitrarily small for L sufficiently large so that condition C2 is satisfied. 

It is easily seen, using DLR structure of the multi-grancanonical Gibbs field, that the case 
when x,y are close from the exterior to the two opposite vertical faces (at distance 3L) can be 
treated exactly like in the homogeneous (constant activity) case; thus we will only consider the 
above mentioned case of x, y belonging to upper and lower sets of conditioning spins. 

Sometimes, just for the sake of simplicity of notation, we will actually drop the explicit 
dependence on the boundary condition r (even though this dependence is crucial). We express H v 
exactly as we did in (3.13): 

H V (V)= E E H B k2 (Pk 2 ) + W Bk2tV \ Bk2 (p k2 \a) 

ki:A kl eAv k 2 :B k2 €l3v 

+ E H c k3 (lk 3 ) +W Ck3 y\ Ck3 ( lk3 \P,a) (j4L7) 
k 3 .c k3 ec v 

+ E H D k4 (h,) + W Dk4 y\ Dk4 (5 k4 \ 7 ,f3,a) 

k 4 :D ki eT> v 

where, as in (3.12), 

W / a 1 ,a 2 (??a 1 |??a 2 ) := W(rj Al \r]A 2 ) = H AlUA2 (r] Al , Va 2 ) ~ H Al (rj Al ) - H A2 (r] A2 ) (A1.8) 

We now proceed to the summation over the 5, 7, (3 variables; we repeat exactly the same operations 
of splitting and gluing that we performed in section 3. We get: 

Z v = E II exp{H(a kl )} [Z Dk± ((0), (a kl ), (0)) Z Dki _ ei ((0), (0), (a kl ))] _1 

a k 1 :A kl €Av 



J] Z Sk2 ((0),(ar,(a) d )J2^(P) 

i k2 eBv 2 p 

n [^((o))]" 1 n (i+<k/3))e^(7) 

"fc 3 GCv k 3 :C k3 eCv 7 

n (i +<(«,/?, 7) ) n 0+*s>)) n [wwr 1 



X 

k4.:D k4 €T>v k4,:D k4 eT>v fc 4 :-Djfc 4 €XV 

(A1.9) 

where the terms [Z Dki _ ei ((0), (a kl ), (0))]" 1 , [Z Dki ((0), (0), (a kl ))]" x (defined in (3.29))j3ome from 
the splitting described in (3.30): in (A1.9), by an abuse of notation, we still denote by C, B and A 
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their truncation in Ri^L- Indeed, since we have generic and not periodic b.c, we have to introduce 
the modifications described in Section 3 (below Proposition 3.2) in fi^P , [ig k as well as in the <3? 
and error terms. Moreover notice that in the expression in (Al.7) above, we continue to denote 
by a, (3, 7, 5 also the configurations on the A, B, C, D blocks outside V; in other words we continue 
to denote by a, /3, 7, 5 also the part of the r (exterior) configuration in A, B, C, D sub-lattices. We 
did not have them in (3.34) since, there, we were using periodic boundary conditions. Now we 
continue with the same operations of splitting as in (3.38) (and gluing as in (3.50) ) only for the 
B (and A) blocks in A, B namely outside the two vertical column V\. It is clear that we cannot 

perform the gluing operation described in (3.38) for the B blocks in V\ and obtain a small value 

(2) (2) 
for the term <J>^ (a). Indeed to get a good upper bound for sup a \&b (a) I we need the validity 

of condition CI, with a sufficiently small e\ for horizontal Rl ,3L rectangles and this condition is 

supposed to hold only for Rl ,3L rectangles completely contained in one of the three squares Q^, 

Q C L or Q r L . For Rl ,3L rectangles centered at B block in B\ we cannot use condition CI since these 

rectangles have simultaneously non-empty overlap with two of the big squares namely Q 1 ^, the 

rectangles Rl„,3L having non-empty overlap with Q c ^, Q 7 ^ behave exactly like in the homogeneous 

case since the activity in Q c ^ U Q 1 ^ is supposed to be constant. 

In this way we obtain the following expression 

a k i: A kl eA 

- n n (i+<<«)) n (i+*k>>) 



ki:A kl eA 



k4,:D k4 eT> v 



P k 3 :C k3 €C v 

xE^(7) II (l + 1>S 4 KA7)) 

7 fc 4 :D fc4 eX>v 

where fJ-A kl (afci) is defined in (3.52), Zy is given by 

~ z v= n z x k m n i^M' 1 n 



(ALIO) 



knA kl eA 



k 2 :B k „eB 



k 3 :C k ^eC v 



n z Dk4 ((o)) 

ki-D k4 eT> v 



and 



Z Vl (a):= J] exp^KJ)^^^),^),^))]" 1 ^^^^),^),^))]- 1 

Ki'-A kl €V t 

a k 2 +e 2 

K 2 :B k GVi 



(ALII) 



Let us call a\ the complex of a variable in A\. If we perform, in the r.h.s. of (ALIO) the sum 
over the 7,/? variables and over the a variables in Ay, we get: 

Z V = Z V Y, ZvMWv&i) (A1.12) 



61 



where, of course, 



(A1.13) 



Sy(a;) = II VA kl (a kl ) 

x n (^^k)) n n o+<>)) 

fci:A fcl e^ k 2 :B k2 €B k 4 :D k4 eV v 

/3 k 3 :C k3 eC v 

xEffW II (l + < 4 («,A7)) 

7 k 4 :D k4 eT> v 

Like in Section 3 we can write: 

n 

Hi R n :RiCV, 1=1 

RjnSj=( ( i<j = l,...,n 

where the polymers Ri are defined like in Section 2 with the obvious changes. In this way we are 
reduced to one-dimensional system on Vi, with finite norm, rapidly decaying interaction. Indeed 
we can write: 

Z r v = exp(tf(a,)) (A1.15) 

Where 

H(a t ) := const. + ^ H(a kl ) - logZJ,^ (a kl ) - logZ^ ^ (a fcl _ ei ) 

+ logZ- (a fc2 ,a fe2+e2 ) + logZ~ fct (Q fer ,4'* ) +e 2 )+ E *r(<*r) 

k 2 --B k2 ev l 2 i rcAi 

where 

1) *f(a r ) := E^,...,^ M^i, ■ ■ ■ , *n) U7=i Cfl>r), 

2) the sum ^ runs over the clusters of (incompatible) polymers "touching" the whole set 
r of ^4-blocks in the sense that the product of the activities of the polymers R±, . . . , R n explicitly 
depend on all the a-variables corresponding to the A-blocks in T and does not depend on any 
other a. 

3) we introduced Z Dk (a k ) = Z Dk ((0), (a k ), (0)), Z~^ (a k , a k+e2 ) = Z~^ ((0), (a k+e2 ), (a k )) 

4) kf is the index of the uppermost 5-block in Vf. 

k* ■= - (M/2 - 1) /2, (M/2 - 1) /2 
anda k T ) +e2 is the configurations in the A-b locks immediately outside (on the top) of V\. 
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Notice that the dependence on the boundary condition r external to V is really present 
(beyond the term Z~ ^ {ctk* , aj^ +e2 )), only in Zo k with Dt adjacent to the boundary dV (upper 

and lower side). 

Prom (A1.15),(A1.16) and the general theory of cluster expansion (see Proposition 3.2) it 
follows that, for M sufficiently large, there exist positive constants ci, mi, m 2 , m 3 such that: 

||^||ooe mi|r| e m2diam r < oo, (A1.17) 

r3A 

and, for any y G d v ,T C Ai,: 

sup sup |$ T (a;) - $ T '(a*)l < cie"" 13 dist ( r 'f) (A1.18) 

a I t,t':t x =t x ,\/ x^y 

We are now reduced to a one-dimensional system with finite norm, rapidly decreasing potential. 
We can then apply the theory developed in [C02] and especially in [CCO] (see also [CEO],[C01]). 

Let us summarize the strategy of [C02], [CCO] to find good mixing properties of the Gibbs 
states for the one-dimensional systems like ours. Consider the system of M/2 variables on 
V\. Suppose that the integers p,n are such that M/2 is a multiple of pn. We divide the interval 
[1, . . . , M/2] into m = intervals ii, . . . ,I m of length pn. We call long range the contribution to 
the interaction coming from the terms with range larger than p. We decompose the potential $ as: 

$ = $ sr + $ lr with $ s r r = if diamT > p; ^ / only if diamT > p (A1.19) 

The idea is to treat & r as a small perturbation. Indeed given a single block Ij, a uniform upper 
bound on the sum of the absolute values of the contributions of the long-range terms involving 
Ij is of the order of nexp(— cp) for a suitable positive constant c. On the other hand for the 
"reduced" system with only short range interactions we can exploit the one-dimensionality and 
the uniform boundedness of the interaction. Indeed the short range transfer matrix has a uniform 
positive gap in its spectrum. This would imply an exponential clustering of the short range Gibbs 
measure: the truncated correlations at the extrema of an interval Ij would decay as exp(— dn) with 
c' depending only on the gap of the transfer matrix. In the perturbative expansions in [C02],[CCO], 
the intervals Ij involved in at least one long range term are treated separately from the other ones 
and they happen to be very rare; on the other ones, where only the short range terms are present 
the mechanism of strictly positive gap of the transfer matrix is active, inducing exponential decay 
of correlations. We refer to [C02], [CCO] for more details; in these articles, (actually in a more 
complicated situation), analyticity of the free energy and decay of truncated correlations are proved. 
In our case, as a consequence of the methods of [C02], [CCO] we get exponential decay of truncated 
correlations. This, together with (Al.18) allows to conclude the proof of Proposition 3.2 □ 

A. 2 A counterexample to USM =>■ MUSM in dimension 3. 

We give here an example that, in general, the implication USM(^l) =>■ MUSM(^l) does not 
hold. We stress that our example is ad hoc, in particular the interaction is translation invariant 
only by even shifts. We believe however it sheds some light on the pathologies that may happen. 
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interacting planes 




e 3 



Fig. 9a 

It is convenient to describe the example (see Fig. 9) by using spin variables, a £ {—1, l} z3 . 
We denote by (ei, e2,es) the canonical basis in Z 3 . The one body potential (magnetic field) is given 
as follows 

-r / \ I — ka x if Xi is even 



ka x if xi is odd 



the two body potential is instead given by 



( —Ja x o y if xi is even, yi is odd and y = x + e\ + ae 2 + bes 
<^ {x , y} (a x ,a y ) = l for some (a, 6) G {(0,0), (1,0), (0,1), (0,-1)} 

I otherwise 

where J > 0. All the other potentials vanish, i.e. <&a = for |A| > 2. 
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Note that the layer {x : x\ = a, a even} interacts only with the layer {x : x\ = a + 1}; in 
particular each double-layer is independent of everything else. Furthermore we claim that each 
double-layer is isomorphic to a standard two dimensional Ising model with staggered magnetic 
field (see Fig. 9b). 




(2x2, x 3 ) if ^3 is even 
(2x2 — 1, X3) if x 3 is odd 



Fig. 9b 

We can in fact map the layer {x : x\ = a, a even} to the even sub-lattice of Z 2 , as follows 

(a,x 2 ,x 3 ) h 

and the layer {x : x\ = a + 1} to the odd sub-lattice of Z 2 , as follows 

, , /(2x 2 -l,x 3 ) if x 3 is even 

(a + l,x 2 ,x 3 ) ^ < } \ - f ■ jj 

[(2x 2 ,x 3 ) if x 3 is odd 

It is easy to verify that under the above mapping the double layer {x : X\ = a, a even} U {x : x\ = 
a + 1} is mapped onto the two dimensional Ising model with the following interaction 



$a(ct) 



— ka x if A = {x} and X\ + X2 is even 

ka x if A = {x} and X\ + X2 is odd 

—Ja x a y if A = {x, x + ei}, A = {x, x + 62} 

otherwise 
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and we are left with studying the strong mixing properties of such a model. 

Let us denote by the Gibbs local specification associated to the interaction (A2.1) and 
by fi^ h the measure obtained from by adding a (constant) magnetic field h €. M., i.e. — k on 
the first line of (A2.1) becames —k — h whereas k on the second line of (A2.1) becames k — h. 
We claim that, if k is chosen large enough (depending on J) such a measure does satisfy condition 
GUSM. Roughly speaking, we have a large magnetic field in either the odd or the even sub-lattice, 
therefore the phase is determined on that sub-lattice; since the other sub-lattice is conditionally 
independent (given the first sub-lattice) we get the strong mixing condition. Indeed one can verify 
that the finite size condition CI holds on squares of side 2 with constants uniform in h. 

On the other hand it is very easy to show that there is no Iq such that the 3 dimensional model 
we started from satisfies GMUSM. Let £ be an odd integer and consider A = Qg {{—{£ — l)/2, 0, 0))U 
Qg ((£ + l)/2, 0, 0)); put a magnetic field hi = —k (resp. /12 = +k) on the first (resp. second) cube. 
The image, under above mapping, of the double-layer {x : x\ = 0} U {x : x\ = 1} is now the 
standard two-dimensional Ising with zero magnetic field. If J is chosen large enough we then have 
a long range order, hence (2.3) fails to hold. 

The pathology that has occurred is the following. Even if the local specification does sat- 
isfy the strong mixing condition separately in each one of the two cubes Qt{{—{£ — l)/2,0, 0)), 
Qt ((£ + l)/2,0, 0)), when we put them together we have a long range order which propagates 
inside the double-layer which sits across the interface between the two cubes. 
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